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Abstract 

In  the  solid  desiccant  wheel  air-conditioning  system,  the  performance  of  the  desiccant  wheel  is 
critical  to  the  capability,  size  and  cost  of  the  whole  system.  Constructing  mathematical  model  is  an 
effective  method  for  analyzing  the  performance  of  the  rotary  wheel  as  well  as  the  system.  The  model 
can  also  be  used  to  guide  system  operation,  interpret  experimental  results  and  assist  in  system  design 
and  optimization.  The  overall  objective  of  this  paper  is  to  provide  a  review  of  various  efforts  that 
researchers  have  made  to  mathematically  model  the  coupled  heat  and  mass  transfer  process 
occurring  within  the  wheel.  The  paper  first  briefly  describes  desiccant  wheel  including  fundamental 
principle,  heat  and  mass  transfer  mechanism  and  the  method  of  model  establishment.  Then  various 
models  consisting  of  ideal  assumptions,  governing  equations,  auxiliary  conditions,  solution  methods 
and  main  results  are  presented.  The  models  can  be  classified  into  two  main  categories:  (1)  gas-side 
resistance  (GSR)  model;  (2)  gas  and  solid-side  resistance  (GSSR)  model  which  can  be  further 
subdivided  into  pseudo-gas-side  (PGS)  model,  gas  and  solid-side  (GSS)  model  and  parabolic 
concentration  profile  (PCP)  model.  It  shows  that  GSSR  models  are  higher  in  precision  and  more 
complex  compared  with  GSR  models.  In  addition,  the  simplified  empirical  models  based  on 
measured  data  are  briefly  discussed.  This  review  is  useful  for  understanding  the  evolution  process  and 
status  quo  of  the  mathematical  model  and  highlighting  the  key  aspects  of  model  improvement  such 
as  taking  account  of  pressure  loss  or  air  leakage. 
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1.  Introduction 

Most  of  traditional  vapor  compression  air  conditioning  system  influences  the 
environment  in  three  aspects.  First,  the  CFCs  and  HCFCs  used  as  working  fluids  are 
ozone-depleting  substances.  Second,  most  of  the  refrigerants  such  as  HCFCs  are  the  global 
warming  gases.  Third,  refrigeration  system  consumes  energy  which  in  general  is  produced 
from  fossil  fuel  and  related  to  release  of  pollutions.  These  polluting  materials  such  as  CO, 
SO*,  NO*  are  harmful  to  the  health  of  people.  While  desiccant  wheel  air  conditioning 
systems  offer  a  lot  of  significant  advantages  over  other  AC  systems.  A  few  remarkable 
advantages  can  be  summarized  as: 

1.  Water,  a  natural  working  fluid,  is  used  as  refrigerant.  Desiccant  materials  such  as  silica- 
gel  and  zeolite  are  also  environmental  friendly. 

2.  Desiccant  cooling  system  can  be  powered  by  low  grade  thermal  energy  sources  such  as 
solar  energy,  geothermic  energy  and  waste  heat.  Hence  the  operating  costs  can  be 
significantly  reduced. 

3.  The  system  deals  with  the  sensible  and  the  latent  loads,  respectively,  and  can  be  utilized 
in  various  climates  or  environmental  conditions. 

4.  It  overcomes  the  discontinuous  problem  of  the  fixed  bed  desiccant  cooling  system. 

5.  Compared  with  the  liquid  desiccant  system  in  which  the  liquid  and  air  directly  interact, 
the  solid  one  is  compact  and  less  subject  to  corrosion. 
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Nomenclature 

a  the  distance  in  radial  direction  (m),  Eq.  (20) 

o' 

A  channel  cross  section  area  (m  ) 

o  o 

Av  transfer  area  per  unit  volume  (m  /m  ),  Eq.  (65) 

B  outside  atmospheric  pressure  (Pa),  Eq.  (26) 

cp  specific  heat  (J/kgK) 

de  hydraulic  diameter  de  =  AA/P,  Eq.  (83) 

De  desiccant  effective  diffusivity  (m2/s),  Eq.  (123) 

Dg  gas  phase  diffusivity  (m2/s),  Eq.  (116) 

Dk  Knudsen  diffusivity  (m2/s),  Eq.  (132) 

Dq  ordinary  diffusivity  (m2/s),  Eq.  (132) 

Z)s  surface  diffusivity  (m2/s),  Eqs.  (108),  (131) 

/  desiccant  fraction  in  the  wall  or  in  the  matrix 

fsptfd  the  mass  of  support  materials/desiccant  per  length  in  a  tube  (kg/m),  Eq.  (32) 

Fs  ratio  of  free  flow  area  to  cross  section  area  of  rotary  wheel,  Eq.  (65) 

h  air-side  convective  heat  transfer  coefficient  (W/m  K) 

hf  compound  convective  heat  transfer  coefficient  (W/m2K),  Eq.  (83) 

hm  convective  mass  transfer  coefficient  (m/s),  Eq.  (47) 

hfg  heat  of  vaporization  of  water  (J/kg) 

H  enthalpy  (J /kg) 

k  thermal  conductivity  (W/mK) 

Ky  gas-side  mass  transfer  coefficient  (kg/m  s) 

Kv  composite  mass  transfer  coefficient  (kg/m2  s),  Eq.  (84) 

L  actual  process  channel  length  (wheel  thickness)  (m) 

Le  Lewis  number 

m  mass  (kg) 

m  mass  flow  rate  of  the  stream  (kg/s) 

m '  rate  of  phase  change  per  unit  length  (kg/m  s),  Eq.  (42) 

r\ 

mf  mass  flow  rate  of  air  per  unit  cross  section  area  (kg/m  s),  Eq.  (65) 
n  the  number  of  the  channel,  Eq.  (20) 

NTU  number  of  transfer  units,  Eq.  (91) 

P  channel  cross  section  perimeter  (m) 
qst  heat  of  sorption  (J/kg  adsorbate) 

r  radial  displace  or  radius  of  the  particle  (m),  Eq.  (143) 

rx  radius  of  the  matrix  axis  (m),  Eq.  (100) 

r2  radius  of  the  rotary  wheel  (m),  Eq.  (100) 

R  separation  factor  (reflecting  desiccant  isotherm  shape),  Eq.  (16) 

T  temperature  (K) 

t  time  (s) 

tr  time  required  for  a  complete  matrix  rotation  (s),  Eq.  (82) 
u  air  velocity  (m/s) 

V  the  total  volume  of  desiccant  wheel  (m  ),  Eq.  (107) 

W  desiccant  adsorption  mass  (kg  adsorbate/kg  adsorbent) 

Xm  actual  channel  width  (m),  Eq.  (20) 

Y  humidity  ratio  (kg/kg) 
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Ya  absolute  humidity  ratio  of  the  air  stream  (kg  water  vapor/kg  dry  air) 

Yd  humidity  ratio  in  equilibrium  with  the  desiccant  (kg  water  vapor/kg  dry  air) 
Yg  absolute  humidity  ratio  of  the  air  stream  (kg  water  vapor/kg  moist  air) 

Ym  humidity  ratio  in  equilibrium  with  the  matrix  (kg  water  vapor/kg  dry  air) 

z  axial  displace  through  matrix  measured  from  period  entrance  (m) 

Greek  letters 

0  relative  humidity 

r\  fraction  of  phase  change  energy  delivered  directly  to  the  air,  Eq.  (42) 

p  density  (kg/m  ) 

o  volume  fraction,  Eq.  (53) 

5  the  thickness  of  desiccant  felt  (m) 

9  time  from  the  beginning  of  period  (s),  Eq.  (82) 

Oj  duration  of  period  j  (s),  Eq.  (82) 
co  wheel  rotate  speed  (r/h) 

t  tortuosity,  Eq.  (109) 

et  felt  porosity,  Eq.  (131) 

Pj  dry  matrix  mass  divided  by  dry  fluid  mass  in  the  matrix  in  period  j,  Eq.  (79) 
Subscript 

0  initial  estate 

1  the  process  air 

2  the  regeneration  air 

a  dry  air 

ad  air  in  desiccant  pore,  Eq.  (116) 

d  desiccant 

eff  effective  property,  Eq.  (140) 

/  entire  felt  (matrix  and  pores),  Eq.  (139) 

fd  dry  solid  felt  material,  Eq.  (138) 
fm  matrix  part  of  felt,  Eq.  (138) 
fp  pores  (inter-particle  voids),  Eq.  (138) 

fs  dry  solid  felt  material,  Eq.  (138) 
g  gas  including  dry  air  and  vapor 

in  inlet  condition,  Eq.  (5) 

j  period  index  (during  the  period  J),  Eq.  (77) 

/  liquid  water  of  frost,  Eq.  (52) 

le  local,  Eq.  (37) 

m  matrix 

s  surface,  in  gas  phase  adjacent  to  gel  particles,  Eq.  (130) 

sp  support  materials,  Eq.  (30) 

v  water  vapor 

vs  saturated  water  vapor,  Eq.  (26) 

w  wall,  Eq.  (119) 

ws  saturated  water,  Eq.  (75) 
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Desiccant  wheel  found  in  the  desiccant  cooling  system  is  widely  acknowledged  as  the  most 
critical  part  of  the  system,  both  in  terms  of  cycle  performance  and  system  cost.  It  has  been 
successfully  used  in  many  systems  for  dehumidification  and  enthalpy  recovery.  With  the 
development  of  computer  hardware  and  numerical  methodology,  mathematical  models  are 
being  used  to  carry  out  critical  investigations  concerning  on  the  desiccant  wheel.  The 
advantages  of  this  method  are  that  it  takes  less  time  and  cost  than  experimental  method  for 
predicting  the  performance  of  a  desiccant  wheel.  The  second  point  is  that  the  mathematical 
model  can  produce  large  volumes  of  results  at  virtually  no  added  expense  and  it  is  very 
convenient  to  perform  parametric  research  and  optimization  analysis.  The  third  reason  is 
that  by  experimental  means,  some  parameters  like  that  inside  the  tubes  cannot  be 
measured.  The  last  but  not  the  least  reason  is  the  fundamental  physics  and  surface 
chemistry  of  rotary  desiccant  wheel,  fixed  desiccant  equipment  and  enthalpy  wheel  are 
similar,  so  that  other  installations  may  be  investigated  too  by  utilizing  the  model  of 
desiccant  wheel  [1-3].  Consequently,  constructing  valid  mathematical  model  for  desiccant 
wheel  has  become  the  key  subject  of  many  studies.  Lots  of  mathematical  models  have  been 
constructed  and  employed  to  analyze,  develop  and  design  desiccant  wheels.  However,  to 
the  authors’  knowledge,  there  is  not  yet  a  work  that  summarizes  comprehensively  the 
numerous  significant  works  done  on  modeling  the  desiccant  wheel  system.  Hence,  the 
objective  of  this  article  is  to  examine  the  progress  in  the  area  of  mathematical  modeling  on 
the  desiccant  wheel. 

Based  on  the  heat  and  mass  conservation  principles,  the  combined  heat  and  mass 
transfer  in  desiccant  wheel  can  be  formulated  with  complex  mathematical  detail  as 
governing  equations  and  auxiliary  relations.  This  is  the  main  and  universal  method  for 
constructing  mathematical  model.  There  are  also  some  simplified  models  for  ease  of 
computation  such  as  pseudo-steady-state  model  [4,5],  which  are  not  fundamentally  correct 
or  universal.  Moreover,  it  is  an  effective  way  to  get  the  empirical  performance  fitted 
correlation  using  large  amounts  of  experimental  results.  But  its  application  is  restricted  in 
the  specific  range  of  experimental  working  conditions.  The  review  of  this  paper  will  focus 
on  the  mainstream  mathematical  models.  Besides,  some  empirical  models  will  be  briefly 
introduced  too. 

The  first  category  of  the  paper  is  distinguished  by  whether  the  model  gives  solid-side 
resistance  (SSR).  Within  the  first  section,  it  is  ordered  according  to  the  dimension  of 
models.  For  the  second  category,  papers  are  grouped  together  by  the  method  of  calculating 
SSR.  For  each  model  reviewed,  attempt  has  been  made  to  gather  up  the  assumptions, 
governing  equations,  auxiliary  relations  and  the  solution  techniques  for  ensuring  the 
integrality  of  the  model.  The  results  and  validation  of  some  models  are  listed  when 
comparison  between  them  are  needed. 

2.  Problem  description 

2.1.  The  fundamental  principle 

Desiccant  wheel  is  also  known  as  dehumidification  wheel,  wheel  dehumidifier  or  rotary 
dehumidifier  in  different  literatures.  The  main  components  of  desiccant  wheel  include 
matrix,  clapboard,  the  wheel  case,  the  air  heater,  the  driven  motor,  etc.  as  shown  in  Fig.  1. 
Matrix  usually  contains  supporting  material  and  desiccant  material.  The  case  of  desiccant 
wheel  is  driven  by  an  engine  and  moves  at  a  given  rotary  velocity.  It  consists  of  a  large 
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Fig.  1.  Schematic  figure  of  rotary  desiccant  wheel. 


number  of  channels  whose  walls  are  constituted  by  supporting  material  and  coated  or 
impregnated  with  desiccant  material.  The  cross  section  of  wheel  is  divided  into  process  air 
side  and  regeneration  air  side  by  the  clapboard.  Water  vapor  is  adsorbed  by  the  desiccant 
when  moist  air  is  passing  through  process  air  side.  And  while  the  hot  air  which  is  heated  up 
by  the  heater  flows  through  regeneration  side,  water  is  desorbed  out  from  the  desiccant  and 
desiccant  material  is  regenerated.  The  airflows  can  be  either  in  counter-flow  or  co-current 
arrangement.  It  can  be  found  that  the  reduction  of  moisture  in  the  process  air  is  at  the  price 
of  the  energy  consumed  to  heat  the  regeneration  air.  This  mode  is  called  active  desiccant 
wheel.  If  no  air  heater  is  built  in  system,  the  mode  is  named  as  passive  desiccant  wheel, 
enthalpy  wheel  or  rotary  energy  wheel.  It  should  be  noted  that  the  wheel  is  usually  installed 
with  thermal  insulation  and  air  proof  material.  Hence  the  energy  and  mass  exchange 
between  the  wheel  and  environment  is  neglected  in  most  cases. 

2.2.  Mechanism  of  heat  and  mass  transfer  in  desiccant  wheel 

During  adsorption,  heat  of  sorption  is  released  as  desiccant  adsorbs  water  vapor.  The 
heat  generated  in  the  desiccant  is  transmitted  through  the  material  which  decreases  the 
sorption  capacity.  Therefore,  the  heat  and  mass  transfer  within  the  desiccant  are  coupled 
and  should  be  considered  simultaneously  in  developing  mathematical  model. 

In  a  desiccant  dehumidifier,  there  exist  gas-side  resistances  (GSRs)  and  solid-side 
resistances  (SSRs).  GSRs  come  from  diffusion  and  heat  conduction  within  the  air  and 
convective  heat  and  mass  transfer  between  air  and  desiccant.  It  influences  the  heat 
and  mass  transfer  from  the  air  stream  to  the  surface  of  the  desiccant  felt.  Due  to  diffusion 
and  heat  conduction  within  the  air  is  small  in  comparison  with  convective  transfer,  they  are 
always  neglected.  SSRs  come  from  heat  conduction  and  mass  diffusion  within  solid 
desiccant  felt.  Where  the  mass  diffusion  involves  a  number  of  physical  mechanisms:  (1) 
surface  reaction  combined  with  the  adsorption  process;  (2)  interparticle  channel  diffusion 
and  adsorption  within  the  layer;  (3)  micropore  and  macropore  diffusion  and  adsorption 
within  the  desiccant  particles  [3]. 
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2.3.  Mathematical  model  establishment 

The  main  steps  to  establish  mathematical  model  include:  (i)  assumption  proposal  and 
selecting  appropriate  control  volume;  (ii)  deriving  the  partial  differential  governing 
equations  based  on  the  mass  and  energy  balances;  (iii)  providing  auxiliary  relation  to  close 
the  governing  equations;  (iv)  adopting  proper  mathematic  approach  to  solve  the  model. 
Some  ideal  assumptions  are  needed  for  the  reason  that: 

(1)  The  mass  and  heat  transfer  occurring  in  the  wheel  are  too  complicated  to  understand 
completely. 

(2)  Some  factors  that  do  not  influence  heat  and  mass  transfer  significantly  can  be 
neglected.  Then  the  complexity  of  the  governing  equations  can  be  reduced  and  the 
solution  is  relatively  simple  and  time-saving. 

2.3.1.  General  governing  equations 

Euler  cylindrical  coordinate  system  has  three  coordinates:  axial  direction  (z),  radial 
direction  (r)  and  circumferential  direction  (cp)  and  is  usually  employed  to  describe  the 
wheel  which  has  the  same  cylindrical  geometry  like  Euler  cylindrical  coordinate  system.  In 
general,  governing  equations  involve  four  main  equations  such  as:  mass  and  energy 
balance  equations  in  the  air  and  in  the  solid  side.  For  clarity,  the  system  is  considered  to  be 
one-dimensional  and  explanations  of  main  terms  in  governing  equations  are  presented 
below.  The  physical  model  is  shown  in  Fig.  2. 

The  moisture  conservation  in  the  air  can  be  expressed  as 

dePa  ^W  +  uEl)=Ky{Yd~Ya)-  { 1 ) 

The  first  term  on  the  left-hand  side  (LHS)  of  the  above  equation  is  the  moisture  storage 
term  in  the  air.  The  second  term  on  the  LHS  represents  the  rate  of  moisture  variation  in  the 
air  due  to  the  axial  flow  of  the  air.  The  first  term  on  the  right-hand  side  (RHS)  of  Eq.  (1) 
expresses  the  rate  of  moisture  variation  in  the  air  caused  by  the  convective  mass  transfer. 
The  energy  conservation  for  the  air  is  described  as 

/ rT~’  rT~’  J  £\2  rT~’  \ 

deCpaPa  [jr)-  +  u-^)  -  -  °  =  h(Td  -  T a)  +  cpvKy(Yd  -  Ya){Td  -  To). 

(2) 


Fig.  2.  Euler  cylindrical  coordinate  system  and  a  model  of  control  volume  of  the  desiccant  wheel. 
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The  first  term  on  the  LHS  of  the  equation  represents  energy  storage  in  the  humid  air.  The 
second  term  on  the  LHS  calculates  the  rate  of  energy  variation  in  the  air  due  to  the  axial 
flow  of  the  air.  The  third  term  on  the  LHS  represents  the  heat  conduction  in  the  air.  The 
first  term  on  RHS  of  Eq.  (2)  expresses  the  convective  heat  transfer  between  the  air  and 
solid  desiccant.  The  second  term  on  RHS  of  Eq.  (2)  indicates  the  sensible  heat  transfer 
between  the  air  and  solid  desiccant. 

The  moisture  conservation  in  desiccant  is  given  by 

/a  w  a2  w\ 

PdS  (~di~  ~  De  EE)  =  Ky(  Yd  ~  Ya) ■  {3) 

The  first  term  on  the  LHS  of  the  above  equation  is  the  moisture  storage  term  inside 
desiccant.  The  second  term  on  the  LHS  shows  mass  diffusion  within  the  solid  desiccant  in 
axial  direction.  The  first  term  on  the  RHS  of  Eq.  (3)  represents  the  convective  mass 
transfer  between  the  air  and  desiccant. 

The  energy  conservation  in  desiccant  can  be  written  as  follows: 

-  (jM) = - Td) + ka  12  ■  Y-)q" 

+  cpvKy(Y  a  —  Y  d)(T  a  —  Td).  (4) 

The  first  term  on  the  LHS  of  the  above  equation  is  the  energy  storage  term  of  desiccant. 
The  second  shows  the  heat  transfer  due  to  heat  conduction  within  desiccant.  The  first  term 
on  the  RHS  of  Eq.  (4)  calculates  the  convective  heat  transfer  between  the  solid  desiccant 
and  air.  The  second  term  on  the  RHS  expresses  the  influence  of  adsorption  heat.  The  third 
term  on  the  RHS  expresses  the  sensible  heat  transfer  between  the  solid  desiccant  and  air. 

It  is  important  to  note  that  in  the  above  governing  equations,  it  is  assumed  that  all  of  the 
adsorption  heat  is  delivered  to  desiccant. 


2.3.2.  Auxiliary  conditions  and  mathematical  methods 

In  order  to  close  the  partial  difference  equations,  some  auxiliary  conditions  have  to  be 
provided.  They  are  boundary  conditions,  initial  conditions  for  transient  analysis, 
adsorption  equilibrium  correlations.  Also,  the  values  of  adsorption  heat  qst ,  convective 
transfer  coefficient  /z,  thermal  conductivity  k  and  moisture  diffusivity  De  in  desiccant 
material  are  need  to  give. 

a.  Boundary  conditions :  Boundary  conditions  describe  the  behavior  of  the  simulation  at 
the  edges  of  the  simulation  region.  The  boundary  conditions,  commonly  encountered  in 
desiccant  wheel  model,  are  usually  in  three  modes: 


(I)  Boundary  conditions  specify  the  values  or  the  normal  derivative  of  parameters  on  a 
surface  [6],  such  as: 

In  the  process/regeneration  air  (i  =  1  for  process  air;  i  =  2  for  regeneration  air) 
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The  insulated  and  impermeable  boundaries  result  in  [7] 
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(7) 

(8) 


(II)  Periodic  boundary  conditions:  The  periodic  boundary  conditions  at  <P  =  0  [8] 

'  Ya(0,z,t)  =  Ya(2n,z,  t), 

Ta(0,z,t)  =  Ta(2n,z,t), 

<  W(0,  z,  t)  =  W(2n,  z,  t),  (9) 

Td(0,z,t )  =  Td(2n,z,t). 

V 

The  periodic  equilibrium  boundary  conditions  expressed  in  the  form  of  limitation 
[9]: 

for  O^z^L 
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(III)  Boundary  conditions  depict  heat  and  mass  transfer  on  a  surface  [10] 

dYn  dW 

PfpDG  (t,  Z,  0)  +  pdDs  (t,  z,  0)  =  0, 


dr 


dr 


(ii) 


dYd  dW 

PfpDG~^r{t,z,ad)  -  pdDs  (t,  z,  ad)  =  Ky[Yd(t,z,ad )  -  Ya{t,z)\, 


dr 


dr 


(12) 


udTdu  „  m  dTw 

kd  ^  (t,  z,  0)  —  pwOCpW 


dr 


dt 


(13) 


dE 

—kd  —r—  (t,  z,  ad)  =  h[Td(t,z,ad)  -  Ta(t,z )]. 
or 


(14) 


b.  Initial  conditions :  Initial  conditions  give  the  value  of  parameters  at  a  given  starting 
time  (t  =  0).  For  example,  the  initial  condition  of  parameter  X  is 

X(t  =  0)  =  X0.  (15) 

X  is  replaced  by  the  detailed  parameters  in  different  models. 
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c.  Equilibrium  adsorption  equation :  The  equilibrium  equation  determines  the  amount  of 
moisture  that  will  be  adsorbed  by  the  desiccant.  Emphasis  should  be  placed  on  the  fact  that 
adsorption  is  a  dynamic  process.  With  different  operation  conditions,  such  as  relative 
humidity  of  the  air  in  equilibrium  with  the  adsorbent  0,  desiccant  temperature  etc., 
desiccant  has  different  adsorption  capacity  and  equilibrium  adsorption  quantity.  Up  to 
now,  through  experiments  to  obtain  the  dynamic  adsorption  data  and  equilibrium 
equation  for  different  desiccants  is  the  main  method.  As  a  result,  there  are  many  types  of 
equilibrium  adsorption  relationships  including:  relative  humidity  of  moist  air  in 
equilibrium  with  the  adsorbent  0  and  separation  factor  R  are  used  to  calculate  the  water 
content  of  the  adsorbent  W\ 

W/Wmax  =  0/[R  +  (1  -  R)0].  (16) 

Some  researchers  consider  that  W  is  the  function  of  0  and  temperature  of  moist  air  in 
equilibrium  with  the  adsorbent  Tj. 

W  =f(Td,  Y cj)  =f'(Td,  <P).  (17) 

Sometimes  it  is  assumed  that  the  W  is  a  function  of  <P  only  and  invariable  with  temperature 
Td: 

W  =/(*)•  (18) 

Besides,  the  following  relationship  is  utilized  in  many  works: 

0=f(W,  Wmax, R, h*, Pws).  (19) 

d.  Mathematical  solution :  Due  to  the  complexity  and  nonlinearity  of  such  coupled  heat 
and  mass  transfer  models,  it  is  difficult  to  obtain  analytical  solutions.  Hence  numerical 
methods  have  been  applied  to  solve  the  partial  differential  governing  equations.  At 
present,  the  most  widely  used  solution  method  for  desiccant  wheel  model  is  finite  difference 
method.  Other  methods  such  as  finite  volume  method  [7]  and  analogy  method  [11,12]  have 
also  been  used  in  some  works. 

3.  Gas-side  resistance  (GSR)  model 

In  early  research,  due  to  the  lack  of  study  on  complex  diffusion  mechanism  inside  solid 
desiccant,  the  heat  and  mass  transfer  processes  from  air  stream  to  desiccant  are  thought  to 
be  controlled  by  GSRs  only.  Gradients  of  temperature  and  concentration  within  solid 
desiccant  are  ignored.  Thus  the  second  terms  on  the  LHS  of  Eqs.  (3)  and  (4)  are  omitted. 
Some  researchers  only  consider  the  transfer  processes  in  axial  direction.  Others  consider 
that  in  axial  as  well  as  circumferential  direction. 

3.1.  One- dimensional  model 

To  save  the  computation  time,  it  is  desired  to  take  into  account  as  less  dimension  as 
possible  in  the  model.  Since  the  variation  in  moisture  as  well  as  temperature  is  more 
significant  in  axial  direction,  it  is  often  chosen  as  main  direction  for  heat  and  mass  transfer 
in  constructing  the  governing  equations.  Usually,  condensation  is  not  considered  in  the 
models.  However,  it  happens  in  certain  conditions.  Hence  some  researchers  include  it  in 
their  models. 
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3.1.1.  Models  without  considering  condensation 

Charoensupaya  and  Worek  [13]  have  derived  a  simple  one-dimensional  dehumidifier 
model  based  on  the  following  ideal  assumptions: 

1.  The  axial  direction  (z)  heat  conduction  and  molecular  diffusion  within  the  desiccant 
coating  were  negligible. 

2.  There  were  no  radial  direction  (r)  temperature  or  moisture  content  gradients  in  the 
matrix. 

3.  The  axial  direction  heat  conduction  and  water  vapor  molecular  diffusion  in  the  air 
were  negligible. 

4.  The  inlet  air  conditions  were  uniform  radially,  but  might  vary  with  time. 

5.  The  tubes  that  make  up  the  rotary  desiccant  wheel  were  identical. 

6.  All  the  channels  were  assumed  to  be  adiabatic  and  impermeable. 

7.  The  effect  of  centrifugal  force  was  neglected. 

8.  The  thermophysical  properties  of  the  dry  air  and  the  properties  of  the  dry  desiccant 
material  as  well  as  the  matrix  were  constant. 

9.  The  pressure  and  velocity  losses  of  the  air  stream  in  the  axial  direction  were  negligible. 

10.  The  heat  and  mass  transfer  coefficients  between  the  air  stream  and  the  desiccant  wall 

were  constant  over  the  entire  dehumidifier. 


It  can  be  inferred  from  the  assumptions  1  and  2  that  this  model  neglects  the  SSRs.  Because 
of  symmetry  characteristic,  the  middle  of  one  tube  could  be  considered  adiabatic  and  half 
size  of  the  channel  was  selected  as  the  differential  control  volume  as  can  be  seen  in  Fig.  3. 

The  channel  was  assumed  to  be  adiabatic  and  impermeable,  hence  the  conservation  of 
moisture  and  energy  between  the  process  air  stream  and  matrix  were 


a  y 
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Conservation  of  moisture  in  the  matrix  was 
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Fig.  3.  A  model  of  dehumidifier  differential  element  for  GSR  model  [13]. 
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where  mm  is  the  total  mass  of  the  matrix  in  the  dehumidifier 


win  i  —  25pmnXwL. 

Conservation  of  energy  in  the  air  was  calculated  by: 


mg_ 

nXm 


I  Wg 

u  dt 


dHA 

0z  ) 


=  2  Ky(Y 


m 


+  2h{T  m 


where  dHm/dW  means  the  differential  heat  of  wetting  and  dHfJ/dYa  is  defined  as 
dHg  _  1  dHm 

dYa  qs,+fdW 


The  subscript  “m”  in  the  above  equations  means  the  matrix  including  nonhygroscopic 
exterior  surface,  support  materials  and  the  desiccant.  Parameters  n ,  Xm  and  L  represent 
holistic  configuration  of  the  wheel.  The  governing  equations  were  solved  subject  to  the 
equilibrium  relation,  Eq.  (16).  It  was  supposed  that  the  variation  of  adsorption  heat  qst  did 
not  influence  system  performance  greatly,  hence  qst  was  taken  as  constant  and  the  average 
value  of  2700kJ/kg  H20  was  used. 

Then  the  models  were  converted  into  nondimensional  form  and  solved  numerically  using 
an  explicit  finite-difference  method.  This  model  had  been  used  to  calculate  an  open-cycle, 
adiabatic,  solid  desiccant  cooling  system  operating  in  the  ventilation  mode.  The  system 
used  an  adiabatic  dehumidifier,  a  regenerative  heat  exchanger  whose  effectiveness  was 
0.95,  and  two  direct  evaporative  coolers  with  effectiveness  of  0.95.  The  effects  of 
dehumidifier  channel  length,  desiccant  mass  fraction  and  desiccant  isotherm  shape  on 
system  performance  in  terms  of  COP  were  investigated.  The  results  showed  that  there  are 
optimum  mass  fraction  and  desiccant  isotherm  shape  which  will  result  in  an  optimum 
COP.  Furthermore,  the  desiccant  cooling  system  operating  in  the  ventilation  mode  can 
yield  COP  as  high  as  1.4. 

Zheng  and  Worek  [14]  have  established  a  transient  one-dimensional  desiccant  wheel 
model  which  has  similar  assumptions  and  control  volume  as  those  of  Charoensupaya  and 
Worek  [13],  with  certain  modifications.  They  omitted  the  storage  terms  dYa/dt  and  dHg/dt 
in  Eqs.  (20),  (21),  (24)  and  adopted  the  enthalpy  of  water  vapor  to  replace  6 Hg/dYa  in 
Eq.  (24). 

The  new  governing  equations  were  subject  to  boundary  condition  of  Eqs.  (5),  (6)  and 
initial  conditions  with  parameter  Yw  and  Tw.  The  relationship  between  humidity  ratio  and 
relative  humidity  was  described  as 


0.62198P,  _  0.62198^ 
B  -  Pv  ~  B/Pvs  -  <Pd  ’ 


where  Pws  presents  the  saturation  pressure  of  water  vapor 


=  exp 


+  ^2  +  hT  +  2,4 T"  T  + 


2i  =  -5800.2206;  22  =  1.3914933;  23  =  -0.04860239, 

A4  =  0.41764768  x  10“4;25  =  -0.14452093  x  10“7;26  =  6.5459673. 
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The  following  empirical  form  was  used  to  give  heat  of  sorption  [15]: 


hfg 


h*(W/Wma)=  lTA/z* 


exp (XW/WmiiX)  -  exp(2) 
1  —  exp(2) 


X  —  constant, 


where  hfg  is  the  heat  of  vaporization  of  water. 
The  equilibrium  equation  was  represented  by 


0  = 


WR 


Wm,x  +  (R-\)W 


Pvs(T) 

Pvs(To) 


ih*- 1 


And  the  enthalpy  of  desiccant  bed  was  expressed  by 


p  iv 

Hm  =  (1  —  f)cpSpTm  +/  cpdTm  +  /  (Hv  —  qst)  d  W 

L  J  o 


The  simulation  aimed  to  study  different  numerical  methods.  They  adopted  an  implicit 
scheme  finite  difference  method  with  variable  step  sizes  in  time  and  space  coordinates  to 
solve  the  governing  equations.  It  has  been  reported  that  compared  with  explicit  scheme, 
the  implicit  scheme  is  unconditionally  stable,  is  of  high-order  accuracy,  and  enables  more 
rapid  calculation  of  dehumidifier  performance  with  much  larger  time  and  space  steps. 
Meanwhile,  the  simulation  program  was  validated  by  two  methods:  (i)  varying  step  sizes  in 
both  space  and  time;  (ii)  comparing  with  other  available  works,  such  as  results  of  Ref.  [13]. 
Both  validation  methods  verified  the  convergence  and  accuracy  of  the  model.  Using  this 
model,  the  effect  of  rotation  speed  on  the  performance  of  the  adiabatic  desiccant  wheel  was 
parametrically  studied,  and  the  optimal  rotational  speed  was  discovered  by  examining 
outlet  humidity  profiles  of  the  adsorption-side  as  well  as  humidity  wave  fronts  inside  the 
desiccant  dehumidifier. 

Moreover  the  same  model  was  used:  (i)  to  optimize  the  wheel  rotational  speed  [16];  (ii)  to 
investigate  the  influences  of  desiccant  sorption  properties,  heat  and  mass  transfer 
characteristic  and  the  size  of  the  wheel  on  the  dehumidification  performance  of  the 
desiccant  wheel  [17];  (iii)  and  to  investigate  the  wheel  performance  at  the  real  operating 
conditions  which  vary  significantly  from  the  design  point  (i.e.,  at  ARI  conditions)  [18]. 

Zhang  et  al.  [6,19,20]  constructed  a  mathematical  model  of  a  honeycombed  desiccant 
wheel.  The  tube  with  honeycombed  geometry  and  differential  control  volume  are  shown  in 
Figs.  4  and  5.  Heat  of  sorption  was  considered  to  be  delivered  to  the  matrix  only.  And  the 
channel  was  assumed  to  be  impermeable  but  not  adiabatic.  Hence,  the  model  consisted  of 
conservation  equations  of  moisture  and  heat  in  the  air  (Eqs.  (31)  and  (33))  and  within 
desiccant  (Eqs.  (32)  and  (34))  separately,  which  are  different  form  the  governing  equations 
(20)  and  (21). 


.  fdYa.  QYA_rr,rv  v\ 
Apa  \dt  +u  dz  )  ~  KyP{  Yd  Ya)’ 


1  3  W 

-fd—  =  KyP(Ya-Yd), 


(32) 
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r 

Air  duct 


Fig.  5.  The  differential  control  volume  in  the  corrugated  duct  [19]. 


/QT  dT  \ 

Apa(cpa  +  Yacpv ) +  u-^j  =  hP(Td  -  Ta)  +  KyPcpv{ 7d  -  Ya)(Td  -  Ta), 

(33) 

1  07^ 

2  d(Cpd  +  WCpl)  +  f  SpCpsp]  ^ 

=  hP(Ta  -  Td )  +  KyP(Ya  -  Y d)qst  +  KyPcpv(Ya  -  Yd\Ta  -  Td).  (34) 

The  adsorption  heat  of  RD  silica  gel  was  calculated  by  the  following  equation  [15]: 

qst  =  hfg[  1.0  +  0.2843  exp(- 10.28  W)\.  (35) 

The  equilibrium  relative  humidity  on  the  surface  of  RD  silica  gel  was  obtained  by  fitting 
fourth-degree  polynomials  to  the  test  data  provided  by  manufacturer  [21] 

0  =  0.0078  -  0.05759  W  +  24. 16554  W2  -  124.478  W3  +  204.226  W4 .  (36) 

The  relationship  between  humidity  ratio  and  relative  humidity  was  described  as  Eq.  (26). 
The  pressure  of  saturated  water  vapor  Pvs  was  given  by  Antonine  equation: 

3816.44 


Pvs  =  exp  23.196 


Td-  46.13 


(37) 
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In  addition,  convective  heat  and  mass  transfer  coefficients  were  determined  by  [22] 


NukkaP 


4  A 


? 


ShkDaP 

V  =  Pa  4A  . 

where 

Nu/c  =  Shjc , 


=  2.302  X  10"5 


P, 


0 


p. 


P0  =  0.98  X  10"5Pfl,  P0  =  256  K, 


where  Da  is  the  mass  diffusion  diffusivity  of  water  vapor  in  the  air. 

The  model  was  solved  by  backward  implicit  difference  scheme  and  validated  by 
experimental  results  of  a  real  honeycombed  paper-based  silica  gel  desiccant  wheel. 
A  good  agreement  was  obtained  between  the  calculation  results  and  the  experimental 
data.  Using  the  model,  the  temperature  and  humidity  ratio  profiles  in  the  wheel  during 
both  the  dehumidification  and  regeneration  processes  were  analyzed  [20].  The  results 
indicated  that  a  hump  curve  of  air  humidity  ratio  along  the  channel  existed  all  the 
time  in  the  regeneration  side.  The  performance  of  desiccant  wheel  was  significantly 
affected  by  the  hump  curve.  In  order  to  improve  the  performance  of  the  desiccant 
wheel,  it  was  essential  to  accelerate  the  hump  moving  from  the  duct  entrance  to  the  duct 
exit  as  quickly  as  possible.  Also  the  effects  of  velocity  and  inlet  temperature  of  regeneration 
air  and  velocity  of  process  air  on  the  hump  moving  speed  were  investigated  too.  It  was 
found  that  increasing  the  velocity  and  inlet  temperature  of  regeneration  air  could 
accelerate  the  hump  moving.  Meanwhile  the  velocity  of  process  air  had  little  influence  on 
the  hump. 

Simonson  and  Besant  [7]  have  presented  a  numerical  model  for  coupled  heat  and 
moisture  transfer  during  adsorption  and  desorption  processes.  Although  proposed  for 
rotary  energy  wheel,  the  model  is  still  listed  here  because:  (1)  desiccant  wheel  and  rotary 
wheel  have  similar  operation  and  function  of  transferring  heat  and  moisture;  (2)  a  new 
factor  rj  is  proposed  in  their  model  to  account  for  the  distribution  of  phase  change  energy 
between  the  desiccant  and  the  air.  The  tube  geometry  cross-section  and  limited  simplifying 
assumptions  are  basically  similar  to  those  of  Zhang  [6]. 

One  tube  is  chosen  as  the  control  volume  (Fig.  6).  The  conservation  of  energy  in  the  air 
is  given  by 

dT  dT 

PgcpgA  -yf  +  uPgCpgA  =  m'qstrj  +  hP(Tm  -  T g),  (42) 

where  rj  is  the  fraction  of  the  phase  change  energy  entering  the  air.  And  then  the 
first  term  on  the  RHS  of  above  equation  calculates  the  heat  of  sorption  delivered  to  the  air 
side. 
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Fig.  6.  Side  view  of  rotary  energy  exchanger  duct  showing  the  aluminum  surface  [7]. 


The  conservation  of  energy  in  the  matrix  (subscript  “m”,  including  aluminum,  desiccant 
and  moisture)  is  described  by 


m 


8  Tjn 

Q)t 


0 

0Z 


al^al 


STm\ 

0Z  ) 


m'qst(l  -  rj)  +  rh'cpi(Tg  -  Tm)  +  hP(T g  -  Tm). 

(43) 


It  is  supposed  that  the  heat  conduction  through  the  matrix  in  the  axial  direction  is 
dominated  by  the  conduction  through  aluminum  surface  layer  (subscript  “al”)  which  is 
expressed  by  the  second  term  on  the  LHS  of  Eq.  (43).  The  first  term  on  the  RHS  calculates 
the  rest  heat  of  sorption  delivered  to  the  matrix.  The  second  on  the  RHS  shows  the  sensible 
heat  transfer  between  the  solid  desiccant  and  air. 

Conservation  of  mass  in  the  air  results  in  two  equations,  one  for  the  dry  air  and  another 
for  water  vapor,  which  is  different  from  the  general  governing  equation  (1): 

A  it +  E  ^ Pi,uA ) =  (44) 


Ei  +  Ez{paU)-{)-  (45) 

Conservation  of  mass  in  the  matrix  is  given  as 

dW 

Pd,dry^d  kP  ?  (46) 

where  “m7”  means  the  rate  of  phase  change  per  unit  length 

m'  =  hmP(pv  -  pv/n),  (47) 

where  pv  and  pv>m  stand  for  water  vapor  density  and  water  vapor  density  on  the  desiccant 
surface.  The  following  equation  is  used  to  calculate  the  value  of  pv>m\ 


P 


v,m 


Pws(Tm)  w 

RvTm  ( Wm-W)C ’ 


0^  W^Wm 


c 

1  +  c’ 


(48) 
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where  Pvs  is  saturation  vapor  pressure  and  calculated  by 


Pvs  =  exp(F), 

'y  +  c2  +  c3t  +  c4t2  +  c5t3  +  c6t4  +  C7  In  T, 
F=i  c 

~^+C9  +  Cl0T+CuT2  +  CnT 3  +  C13  111  T, 

<  i 


173  K<T <273  K, 
273  K<  F <473  K. 


Heat  of  sorption  is  given  by 

JT  =  h*{ W/  W max)  =  1  +  Ah* ( 1  -  ,  X  =/( IF),  (50) 

hfg  \  WmaJ 

where  X  describes  the  shape  of  heat  of  sorption  as  a  function  of  moisture  content. 

The  geometry  and  property  relations  for  the  moist  air  ( g ),  desiccant  ( d)  and  matrix  ( m ) 
are 


CpgPg  —  Pacpa  T  PvCpv> 

(51) 

Pd  Pd, dry(^  T  1), 

cpdPd  ~  Pd,dry(WCpl  +  Cp£ ry), 

(52) 

Pm  ® dPd  "E  ^alPal’ 

Cpm Pm  =  Pm  =  ®dPdcpd  T  Oa\Pa\Cpal-> 

(53) 

+  °al  —  1  • 

(54) 

The  convective  heat  transfer  coefficient  h  in  Eqs.  (42),  (43)  was  obtained  from  numerical 
simulations  of  simultaneously  (thermal  and  hydrodynamic)  developing  flow  in  an 
equilateral  triangular  duct  [23].  Convective  mass  transfer  coefficient  hm  in  Eq.  (47)  was 
determined  by  using  the  analogy  between  heat  and  mass  transfer.  The  boundary  conditions 
for  the  problem  included  Eqs.  (5)  and  (7)  and 

Pa,m  —  P/,in?  Pv, in  ~  Pi,v, in?  Uw  =  Ui, in-  (55) 

The  governing  equations  were  discretized  by  finite  volume  method  with  a  staggered  grid 
[24].  Upwind  difference  scheme  and  central  differencing  were  used  in  the  air  and  matrix, 
respectively.  Transient  terms  and  source  terms  were  treated  in  fully  implicit  and  semi- 
implicit  manner.  The  discretized  equations  were  solved  using  a  Gauss-Seidel  iteration 
technique  with  under  relaxation. 

The  model  was  validated  by  comparing  its  results  with  the  experimental  data  of  a 
commercial  molecular  sieve  rotary  energy  exchanger  by  Simonson  and  Besant  [25].  It  has 
been  reported  that  the  measured  and  simulated  results  agreed  well  within  experimental 
uncertainty.  Besides,  the  effects  of  some  important  factors,  such  as  phase  change  term  r\ , 
adsorption  heat  and  isotherm,  axial  heat  conduction  in  matrix,  storage  in  air  and  entry 
length,  on  the  predicted  performance  of  molecular  sieve  energy  wheel  were  discussed.  The 
predicted  effectiveness  was  found  to  be  greatly  influenced  by  the  value  of  r\. 

3.1.2.  Models  considering  condensation 

In  the  above  models,  the  desiccant  wheels  are  formulated  at  or  near  atmospheric 
pressure,  where  no  condensation  happens.  When  desiccant  wheel  works  under  high 
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pressure,  condensation  might  occur  in  the  regeneration  process.  The  models  reviewed  in 
this  part  are  capable  of  properly  simulating  desiccant  wheel  operating  at  high  pressures. 

In  order  to  analyze  the  performance  of  an  adiabatic  desiccant  wheel  under  high  process 
and  regeneration  air  stream  pressures,  Mihajlo  and  William  [26]  have  derived  a  model 
containing  two  groups  of  governing  equations,  one  for  the  noncondensation  situation  and 
another  for  the  case  when  condensation  occurred.  The  ideal  assumptions,  governing 
equations  for  noncondensation  situation  and  auxiliary  relations  are  similar  to  that  of 
Zheng  and  Worek  [14].  A  new  set  of  governing  equations  were  constructed  when 
condensation  happened  such  as: 

Conservation  of  moisture  in  the  process  air. 

/  v  v  \  2h{T  a  —  T m) 

PaapU  ^  —  z.Ky{Y m  Y a)  ^  .  (56) 


Conservation  of  moisture  in  the  matrix : 


fPm?> 


8W 

dt 


=  2  Ky{Y  a 


v  \  i  Tm) 

H  m)  i  t 

hfg 


The  temperature  of  air  which  was  saturated  with  water  vapor  was  taken  as  constant 


Ta  =  constant. 


Conservation  of  energy  in  the  matrix 


2tfstKy(  Ya 


v  \  i  2h(T a  Tm) 
Ym)  +  qst - 7 - 

hfg 


An  implicit  finite-difference  scheme  was  adopted  to  solve  the  model.  The  performance 
analysis  using  this  model  revealed  that:  (1)  the  value  of  optimal  separation  factor  R  under 
high  pressure  is  greater  than  that  under  atmospheric  pressure;  (2)  if  the  regeneration 
temperature  Treg  is  between  65  and  100  °C  and  R  is  between  0.01  and  0.1,  condensation 
could  occupy  21-40%  of  regeneration  section  of  the  wheel;  (3)  only  for  R  larger  than  0.4 
does  condensation  not  happen. 

A  similar  extension  of  the  model  [6]  was  conducted  by  Zhong  [27],  to  include 
condensation  process.  Similarly,  two  sets  of  governing  equations  were  established, 
respectively,  for  the  condensation  process  and  noncondensation  process.  The  assumptions, 
governing  equations  of  noncondensation  status,  auxiliary  relations  and  solution  method 
were  similar  to  that  of  Zhang  [6].  The  derived  model  of  condensation  process  consisted  of: 

Conservation  of  moisture  in  the  process  air. 


T  u 


=  KyP(Yd 


Ya)  + 


hP(Td  -  Ta) 

hfg 


Conservation  of  moisture  in  the  desiccant'. 


Yd)  + 


hP(Ta  -  Td) 

hfg 


(61) 
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Conservation  of  energy  in  the  process  air : 


Apa(cpci  “I-  YaCpv) 


dT 


a 


pv , 


dt 


H-  U 


dT 


a 


dz 


=  hP(Td  -  Ta)  +  KyP cpv( Y d  -  Ya)(Td  -  Ta). 


Conservation  of  energy  in  the  desiccant'. 


1  07" 

2  d(Cpd  T  We pl)  +  f  SpCpSp]  ^ 

=  hP(Ta  -  Td)  +  KyP(Ya  -  Y g)qst 

N  hP(Ta-Td ) 

+  KyPCpfYa  -  Y  d){T  a  -  Td)  +  - -fe. 


h 


lfg 


This  model  was  carried  out  to  discuss  the  performance  of  RD  silica  gel  desiccant  wheel 
under  high  pressure.  It  has  been  reported  that  when  the  separation  factor  R  is  0.01, 
operating  pressure  almost  had  no  effect  on  desorption  capacity,  but  sorption  capacity 
enhanced  with  increasing  pressure.  Especially  if  the  pressure  is  higher  than  3  atm,  the 
sorption  capacity  could  be  greatly  improved.  When  R  reaches  1,  the  adsorption  capacity  is 
enhanced  and  desorption  capacity  gets  worse  with  pressure  rising.  If  the  pressure  is  higher 
than  3  atm,  the  variation  of  adsorption  capacity  is  not  significant.  It  was  also  found  that 
desiccant  wheel  achieves  a  better  performance  at  the  condition  of  R  =  1  as  compared  to 
R  =  0.01. 

Condensation  and  frosting  processes  occurring  in  energy  wheels  have  been  investigated 
too.  Simonson  and  Besant  [28]  expanded  the  model  proposed  by  them  [7]  to  include 
condensation  and  frosting  processes.  The  difference  between  the  governing  equations  of 
these  two  models  is  that  the  adsorption  heat  qst  in  Eqs.  (42)  and  (43)  is  replaced  by  heat  of 
vaporization  hfg.  In  saturation  condition,  the  density  of  water  vapor  can  be  approximated 
by 


P 


Pv  P  VS 


vs 


vs  n  'r 

g 


where  Rv  is  the  specific  gas  constant  of  water  vapor. 

The  numerical  schemes  are  similar  to  those  used  in  Ref.  [7].  Mathematically  and 
physically  the  problem  becomes  more  complex  when  condensation  and  frosting  are 
included.  Hence  these  conditions  require  a  different  numerical  solution  strategy  and  then 
different  algorithms  are  adopted  during  sorption  and  saturation  process  as  shown  in 
Table  1.  This  model  could  simulate  the  condensation  and  frosting  happening  in  energy 
wheel.  The  solutions  were  validated  by  comparing  the  simulated  results  with  experimental 
data  on  a  commercial  energy  wheel.  This  study  analyzed  the  sensitivity  of  condensation 
and  frosting  to  wheel  speed  and  desiccant  type.  The  simulation  results  showed  that  a 
desiccant  with  a  linear  sorption  curve  has  smaller  amounts  of  condensation/frosting 
compared  with  Type  I  sorption  curve. 

Dimensionless  equations  of  air-to-air  energy  wheel  were  derived  using  this  model  [29].  It 
was  indicated  that  the  performance  of  energy  wheel  is  a  function  of  the  operating 
temperature  and  humidity.  In  addition,  based  on  this  model,  Jeong  and  Stanley  [30] 
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Table  1 

Different  algorithms  used  to  solve  the  governing  equations  [28] 


Sorption  ($<  100%) 

Common 

Saturation  (E  =  100%) 

1 .  Estimate  the  rate  of  phase  change  and 
the  properties  needed  in  the  governing 
equations 

2.  Solve  the  Tm,  pv,  W,  u  and  TCJ 

2.  Solve  the  Tm,  mr,  W ,  u  and  Tg 

fields  in  order  with  Eqs. 

fields  in  order  with  Eqs. 

(42)-(46) 

3.  Update  the  properties 

(42)-(46) 

4.  Update  the  rate  of  phase 

4.  Update  the  rate  of  phase 

change  (rh!) 

5.  Return  to  step  2  and  iterate  until  a 
converged  solution  is  reached 

6.  Adjust,  as  necessary,  the  algorithm 
used  in  each  grid  for  the  next  time-step 

7.  Increment  time,  return  to  step  2  and 
iterate  until  a  quasi-steady  solution  is 
obtained 

change  ( pv ) 

Table  2 

Features  of  gas  side  resistance  (GSR)  model  reviewed  in  present  work 


Ref 

Dimension 

Condensation 

Boundary  conditions 

Mode  Expressions 

Parameters  of  initial  condisions 

Other  relations 

[13] 

I 

— 

— 

— 

Eqs.  (16),  (25) 

[14] 

I 

I 

Eqs.  (5),  (6) 

T 

Yw 

Eqs.  (26)-(30) 

[6] 

I 

I 

Eqs.  (5),  (6) 

T 

±  a 

Ya  Td  Yd  W 

Eqs.  (26),  (35)— (4 1) 

[7] 

I 

Eqs.  (5),  (7),  (55) 

— 

Eqs.  (47)-(54) 

[26] 

I 

V 

I 

Eqs.  (5),  (6) 

T 

Yw 

Eqs.  (26)-(30) 

[27] 

I 

V 

I 

Eqs.  (5),  (6) 

T 

J  a 

Ya  Td  Yd  W 

Eqs.  (26),  (35)— (41) 

[28] 

I 

V 

— 

Eqs.  (5),  (7),  (55) 

— 

Eqs.  (47)-(54),  (64) 

[8] 

II 

I 

Eqs.  (5),  (6) 

T 

J  a 

Ya  Td  W 

— 

II 

Eq.  (9) 

developed  practical  enthalpy  wheel  effectiveness  correlations  which  can  be  more  easily 
used  to  design  and  analyze  passive  wheel. 

3.2.  Two-dimensional  model 

Multi-dimensional  mathematical  models  reflecting  the  effects  in  radial  or  circumferential 
directions  could  be  used  to  analyze  heat  and  mass  transfer  processes  comprehensively  with 
improved  accuracy.  But  the  complexity  increases  simultaneously.  Usually  two-dimensional 
model  is  recognized  as  a  compromise  method  between  accuracy  and  complexity. 

Yu  et  al.  [8]  have  derived  a  two-dimensional  GSR  model.  Convection  items  dY /dcp , 
dW /dcp  and  dT /dcp  in  circumferential  direction  which  reflect  the  heat  and  mass  transfer 
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due  to  the  rotation  of  the  wheel  were  added  into  the  model.  The  control  volume  of  the 
model  was  the  same  as  that  shown  in  Fig.  2.  The  governing  equations  were  as  follows: 
Conservation  of  moisture  in  the  air. 


87  97  m;.dYa  KvAv 

8 1  8 <p  +  PiFs  9z  PiFs  1  J 

Conservation  of  moisture  in  the  desiccant. 


9  W  9  W 
T  co 


K 


y 


(  7(,  —  7 (/). 


9 1  9  <P  Pcfi 

Conservation  of  energy  in  the  air  equation'. 

9 Ta  9 Ta  m'f  9 Ta  hAv 

+  oo— - b 


9 1  9 cp  PiFs  dz  PjF fcpa  +  Yacpv ) 

Conservation  of  energy  in  the  desiccant. 


(Td  -  Ta). 


9  Td  9  Td 

^r+“— 


h 


9  q>  Pdd(cpd  +  Wcpi) 


(Ta  -  Td)  + 


Kv(Ya  -  Y d)q 


St 


Pdb(cpd  +  Wcp! ) 


Boundary  conditions  of  modes  I  and  II  were  introduced  to  solve  the  problem.  The 
governing  equations  were  discretized  into  finite  difference  equation  group.  For  space 
varying  items,  one  order  upwind  difference  scheme  was  adopted  for  convection  terms.  For 
time  variation  items,  different  schemes  were  used  according  to  the  demands  of 
computation.  To  compute  steady  state,  fully  implicit  scheme  was  employed  because  of 
its  reliability  of  convergence.  To  investigate  the  dynamic  behavior  of  the  dehumidification 
device,  explicit  scheme  which  does  not  need  the  iteration  process  was  recommended. 
Owing  to  the  constraints  of  stability,  time  step  should  be  small  enough,  or  the  divergence 
of  solution  cannot  be  avoided.  Moreover,  small  time  step  was  also  numerically  necessary 
for  accurate  computation.  Crank-Nicolson  scheme  could  be  carried  out  for  dynamic 
investigation  too.  A  proper  mesh  size  should  be  specified  to  save  computation  time.  It  was 
found  that  the  predicted  result  of  this  model  is  more  close  to  the  actual  process  than  that  of 
Charoensupaya  and  Worek  [13].  This  is  mainly  because  although  the  convection  items  in 
circumferential  direction  can  be  neglected  under  low  rotation  speed,  they  influence  heat 
and  mass  transfer  process  in  almost  all  occasions.  Hence,  when  they  are  added  to  the 
model,  more  practical  results  can  be  achieved  compared  with  one-dimensional  model.  Also 
the  analytic  results  indicated  that  when  the  density  of  the  adsorbent  p^is  smaller  than  0.2, 
the  influences  of  dY /dtp  and  dT /dcp  cannot  be  negligible. 

Feng  et  al.  [31]  applied  analogy  method  to  solve  the  proposed  governing  equations  by 
Yu  et  al.  [8].  Their  study  have  shown  that  the  performance  of  dehumidifier  is  a  function  of 
the  nondimensional  structure  parameters  and  thermophysical  property  parameters  (such 
as  nondimensional  rotation  period  t'r  =  (ma*  tr)/(pa*  Fs*  L),  nondimensional  moisture 
ratio  m'  =  ( pa*Fs  •  Ya)/(  Wm^x  •  Mw))  with  the  operation  parameters  keeping  constant.  The 
GSR  models  summarized  in  this  part  are  reviewed  in  Table  2. 


4.  Gas  and  solid-side  resistance  (GSSR)  model 

The  previously  mentioned  GSR  models  do  not  consider  the  heat  conduction  and  mass 
diffusion  within  solid  side,  thus  they  do  not  reflect  the  actual  transfer  process  occurring  in 


o 

ON 


Table  3 

Features  of  gas  and  solid  side  resistance  (GSSR)  models  reviewed  in  present  work 


Refs. 

Dimension 

Conduction  in 

Diffusion  in 

Boundary 

Parameters  of 

Other  relations 

solid-side 

solid-side 

conditions 

initial  conditions 

z 

r 

(p 

z 

r 

(p 

Mode 

Expressions 

[32] 

I 

PGS  model 

I 

Eqs.  (5),  (6) 

— 

Eqs.  (26),  (75),  (76) 

II 

Eq.  (9) 

[9] 

I 

PGS  model 

— 

— 

— 

Eqs.  (28),  (83)— (88) 

[33] 

I 

PGS  model 

0-NTU  Model) 

II 

Eq.  (10) 

YaHg 

Eqs.  (93)— (95) 

[34] 

I 

PGS  model 

I 

Eqs.  (5),  (6) 

T  a  Ya 

Eqs.  (35),  (96),  (98),  (105),  (106) 

(Radial-flow  wheel) 

II 

Eqs.  (103),  (104) 

[35] 

I 

V 

V 

— 

— 

— 

Eqs.  (35),  (96),  (109),  (112)-(114) 

[10] 

II 

V 

V 

Ill 

Eqs.  (1 1 )— ( 1 4) 

Td  Yd 

Eqs.  (11 9)— ( 121) 

[36] 

II 

V 

V 

V 

V 

I 

Eqs.  (5),  (6) 

T a  Ya  Td  W 

Eqs.  (26),  (36),  (126)-(128) 

II 

Eq.  (9) 

[37] 

II 

V 

V 

V 

V 

I 

Eqs.  (5M8) 

— 

Eqs.  (129),  (132),  (135)-(137) 

III 

Eqs.  (12),  (14) 

[38] 

I/II/III 

V 

V 

V 

V 

V 

V 

I 

Eqs.  (7),  (8) 

— 

Eqs.  ( 140)— ( 142) 

III 

Eqs.  (12),  (143)— (145) 

[39] 

I 

PCP  model 

I 

Eqs.  (153),  (156) 

W  Ts 

Eqs.  (36),  (126)— (128),  (158)-(161) 

(Diffusion  in  the  radial 
direction  of  a  particle) 

III 

Eqs.  (154) 

[40] 

I 

PCP  model 

II 

Eqs.  (166),  (167) 

WHgHd  Ya 

Eqs.  (162),  (164),  (165) 
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the  desiccant  wheel.  To  solve  this  problem,  models  considering  heat  and  mass  transfer  in 
both  solid  and  air  side  are  proposed.  Such  models  are  called  GSSR  models  (listed  in 
Table  3).  GSSR  model  is  more  preferable  compared  with  GSR  model  for  it  can  explain  the 
effect  of  heat  and  mass  diffusion  within  solid  desiccant  and  the  precision  of  governing 
equations  describing  the  physics  of  desiccant  wheel  is  enhanced. 

Typical  diffusion  theory  holds  that  water  vapor  molecule  has  three  types  of  diffusion 
mode  in  porous  mediums  which  includes:  ordinary  diffusion,  Knudson  diffusion  and 
surface  diffusion.  The  actual  diffusion  process  is  the  combination  of  these  three 
mechanisms. 

Ordinary  diffusion  which  is  also  called  molecular  diffusion  or  Fickian  diffusion 
occurs  when  the  mean  free  path  is  relatively  short  compared  to  the  pore  size.  This 
diffusion  mode  is  applicable  to  Brownian  motion,  where  the  movement  of  each 
particle  is  random  and  not  dependent  on  its  previous  motion.  The  transfer  diffusivity 
D0  relates  the  macroscopic  flux  of  molecules  in  a  system  to  a  driving  force  in  the 
concentration.  Knudson  diffusion  occurs  when  the  mean  free  path  of  gas  molecules 
is  relatively  long  compared  to  the  pore  size.  It  describes  situations  in  which  gas  mole¬ 
cules  collide  more  frequently  with  flow  boundaries  than  with  other  gas  molecules. 
Knudsen  diffusion  is  dominant  for  pores  that  range  in  diameter  between  2  and  50  nm. 
Surface  diffusion  is  also  used  to  explain  a  type  of  pore  diffusion  in  which  solutes 
adsorb  on  the  surface  of  the  pore  and  hop  from  one  site  to  another  through  inte¬ 
ractions  between  the  surface  and  molecules.  All  of  these  transport  processes  can  be 
described  by  means  of  Fick’s  law  of  diffusion  with  appropriate  diffusion  coefficients 
(D0,  Dk,  T>s). 

Pesaran  and  Mills  [39]  investigated  the  diffusion  mechanism  in  silica  gel.  It  has  been 
reported  that  surface  diffusion  is  the  dominant  one  for  moisture  transfer  in  micro  porous 
silica  gel  such  as  RD  gel  (regular  density  silica  gel).  While  for  macro  porous  silica  gel,  for 
example  ID  gel  (intermediate  density  silica  gel),  both  Knudsen  and  surface  diffusions  are 
important  mechanism.  Ordinary  diffusion  can  be  negligible  for  the  silica  gel  at  atmospheric 
pressure. 

Based  on  different  methods  to  account  for  the  SSR,  GSSR  models  are  classified  into 
three  subcategories:  (i)  pseudo-gas-side  (PGS)  models  in  which  lumped  convective 
coefficients  are  adopted  to  account  for  SSR;  (ii)  gas  and  solid-side  (GSS)  models  in 
which  the  second  order  diffusion  terms  are  added  to  account  for  SSR;  (iii)  parabolic 
concentration  profile  (PCP)  models  which  assume  that  a  parabolic  concentration  profile 
for  moist  concentration  exists  at  all  the  times  in  the  desiccant  particle. 


4.1.  Pseudo-gas-side  (PGS)  model 

A  simplified  method  to  include  the  heat  and  mass  transfer  in  the  solid  desiccant  is  using 
lumped  transfer  coefficients  ( h '  as  heat  transfer  coefficient  and  Ky  as  mass  transfer 
coefficient)  which  include  heat  and  mass  transfer  resistances  not  only  in  gas  side  but  also  in 
solid  side.  Such  modified  models  are  named  as  PGS  model.  Since  the  second-order 
diffusion  terms  are  not  added,  the  complexity  and  computational  time  of  PGS  model  do 
not  increase  compared  with  GSR  model.  The  lumped  transfer  coefficients  are  usually 
obtained  by  the  correlations  of  experimental  data.  The  correlations  for  silica  gel 
as  adsorbent  given  by  the  Hougen  and  Marshall  [41]  based  on  Ahlberg’s  experimental 
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data  [42]  are  in  wide  use: 

h'  =  0.683m;  Re"a51cpa, 

(69) 

r  =  0.704m;  Re-0'51 . 

y  1 

(70) 

4.1.1.  Axial-flow  PGS  model 

Maclaine-Cross  [43]  has  presented  a  finite-difference  model  known  as  MOSHMX 
(method  of  solving  heat  and  mass  exchangers).  The  gas-side  controlled,  lumped 
capacitance  mass  transfer  coefficient  was  used  in  the  model  which  is  actually  a  typical 
PGS  model.  A  model  of  desiccant  wheel  with  PGS  coefficients  similar  to  the  one  proposed 
by  Maclaine-Cross  [43]  has  been  described  by  Holmberg  [32].  In  the  model,  air  thermal  and 
moisture  storage  capacities  dY/dt  and  dT /dt  were  neglected  due  to  their  small  value 
compared  to  the  matrix  capacities.  The  control  volume  is  depicted  in  Fig.  7. 

Mass  and  energy  conservation  equations  were : 


QYg 

dz 


+ 


flnmdW 
L  dt 


dH 


m 


a 


a 


dz 


l  iRm  ^Hm  r  A  d  Tm 


=  0. 


Mass  and  energy  transfer  equations  were: 
.  dYa 


m 


a 


dz 


K'mm 

_ (Y  —  Y  ) 

m  -1  a) : 


dH. 


m 


a 


a 


dz 


K\4 

L 


h 


K 


(Tm  —  Tf)  +  hfg  ( Y m  —  Y  a) 


y 


(74) 


where  it  is  assumed  that  the  rotor  matrix  is  composed  of  supporters  (, sp )  and  adsorbent  ( d)  with 
the  mass  fractions  1  —f  and  f  respectively.  Thus  the  enthalpy  of  the  matrix  is  defined  by 

Hm  =  (1  -f)Hsp  +f[Hd  +  W(Hws  -  qst  +  hfg)\.  (75) 

The  equilibrium  isotherm  was  in  the  form  of  Eq.  (18) 

W  =  0.550.  (76) 


Fig.  7.  The  model  of  control  volume  and  a  side  view  of  one  of  ducts  for  axial-flow  PGS  model  [32], 
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Holmberg  adopted  implicit  Crank-Nicholson  scheme  for  the  energy  conservation 
equation  (72)  and  ordinary  central  difference  schemes  for  other  equations.  Gauss-Seidel 
solution  technique  was  employed  to  solve  the  model.  A  parametric  study  was  conducted  to 
show  the  influences  of  matrix  heat  and  moisture  capacities  on  temperature  and  humidity 
efficiencies.  It  was  indicated  that  there  is  an  optimum  drying  humidity  efficiency  with 
variable  heat  capacities  and  the  optimum  value  increases  in  magnitude  with  increasing 
moisture  capacity. 

Based  on  the  studies  of  Maclaine-Cross  [43]  and  Holmberg  [32],  Jurinak  and  Mitchell  [9] 
neglected  the  axial  heat  conduction  inside  the  matrix  and  developed  a  dimensionless  PGS 
model  in  period  j  (j  is  a  period  index)  of  rotary  dehumidifier.  The  control  volume  was  the 
same  as  shown  in  Fig.  7  and  the  governing  equations  were  given  as: 

Mass  and  energy  conservation  in  period  j : 


QYa 

dz' 

PE 

dz' 


(77) 

(78) 


where  coefficients  pj  and  Kj  are  defined  as: 


Pj  —  ^dry  matrix  mass,  y'/^dry  air  containned  in  matrix,  j 


Mass  and  heat  transfer  equations  in  period  j : 


QYg 
dz ' 


dH, 


—  Awj(Le0cpq(Tm  —  T  a)  +  Hv(Ym  —  Ya )). 


where  dimensionless  coordinates  in  the  governing  equations  are 
zf  =  z/L,  6'  =  9/tr,  0j  —  9j/tr 
and  dimensionless  heat  and  mass  transfer  length  are  given  by 

4ti 


Ajrj  =  JnL/ Vi,  J r  i  — 


Tj 


j> 


TJ 


PadeCpg 


(81) 

(82) 

(83) 


4  K! 

y 

pade 


Overall  Lewis  number  was  expressed  as 


Le  o  —  J  rj/Jwj- 


The  enthalpy  of  the  matrix  at  any  temperature  and  water  content  with  respect  to  a  dry 
matrix  at  273  K  were  given  by 


Hm  —  cpm(Tm  —  273)  +  cpi  W(Tm 


273)  +  hfg 


w 

(1  —  h*)  dW, 


dst_ 

hfg 


h*(W/Wmax)=  lTA/z* 


exp(2IF/IFmax)  -  exp (2) 
1  —  exp(2) 


2  =  constant, 
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where  the  dry  matrix  thermal  capacities,  cpm,  is  related  to  the  specific  heat  of  desiccant  cpd 
and  of  the  nonsorbing  matrix  materials  cpi 

Cpm  =  cpd  +  ^  JjnjCpj) f f .  (87) 


Six  different  isotherm  equations  were  adopted  in  the  study  to  investigate  their  influences. 
They  are  functions  of  W/Wmax  and  here  only  the  simplest  linear  style  is  listed  in  detail: 


&=G(W/  WmSLX)  =  W/  Wmax.  (88) 

A  finite  difference  solution  was  used  to  study  the  effects  of  different  matrix  properties  on 
the  performance  of  a  rotary  dehumidifier,  and  in  each  case  the  results  were  obtained  for  the 
maximum  dehumidification  efficiency.  The  results  were  useful  to  select  desiccant  and 
analyze  the  open-cycle  systems. 

Cejudo  et  al.  [44]  validated  the  model  proposed  by  Jurinak  and  Mitchell  [9]  by 
experimental  data.  Moreover,  they  used  the  real  data  to  train  a  neural  network  to  calculate 
the  parameters  of  a  new  neural  network  model.  Their  results  demonstrated  that  the  new 
neural  network  model  is  more  adequate  to  simulate  wheel  performance  when  heat  losses  of 
the  rotating  wheel  are  taken  into  account. 

Later,  Van  den  Bulck  et  al.  [33]  pointed  out  that,  the  model  solution  in  Ref.  [9]  involved 
considerable  computation  and  was  not  suitable  for  long-term  simulation  of  desiccant 
cooling  systems.  In  order  to  reduce  the  calculation  efforts,  they  introduced  number  of 
transfer  units  (NTU)  in  the  governing  equations.  The  new  conservation  equations  using 
NTU  for  period  j  were 


dYa  dW 

— -  H - 

dz'  dt' 


dHg  dHm 
dz'  dt' 


=  NTU1,J(7m  -  Yd 


(91) 


dM  dH 

9  =  NTU -  Ta)  +  NTU WJHv(Ym  -  Ya), 


dz' 


07 


(92) 


9 


where  air  stream  mass  flow  rates,  wheel  rotation  speed  and  mass  of  desiccant  are  combined 
into  two  dimensionless  capacitance  rate  parameters  r f. 


r,= 


Kid 

trrhj 


(J  =  1,2),  t’  =  il 


0j  r, 


Ossa's; 


l 

h 


z 


'  =  z/L 


(93) 


and  the  number  of  transfer  units  is  defined  as 


NTU.ty  = 


K'A 


m 


(94) 


7 


NTU  tJ  = 


h'Aj 

CRjCpg 


dHd 

ZTjy' 
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This  approach  is  called  effectiveness-number  of  transfer  units  (e-NTU)  method  which 
clearly  represents  the  influence  of  flow  unbalance  and  transfer  coefficients  respectively. 
Boundary  conditions  in  mode  II  and  thermophysical  property  parameters  obtained  from 
Ref.  [45]  were  applied  to  close  the  partial  difference  equations.  Then  the  model  was  solved 
by  means  of  characteristics  and  the  shock  wave  method.  Using  this  model,  the 
performance  correlation  of  ideal  dehumidifier  with  infinite  heat  and  mass  transfer 
coefficients  was  calculated  as  a  function  of  NTU  and  operating  parameters  Ty. 

Subsequently,  correlations  for  the  humidity  and  enthalpy  effectiveness  of  actual  silica  gel 
dehumidifiers  with  finite  transfer  coefficients  as  functions  of  NTU  were  conducted  by  Van 
den  Bulck  et  al.  [46].  In  their  analysis,  equilibrium  isotherm  and  enthalpies  of  gas  and 
matrix  were  calculated  by 

<P  =  (2.\\2W)qJhf»(29.9\Pvs)qJhf«-\  (96) 

Hg  =  cPg(Tg  -  7W)  +  Hv[Tn f]  Yg,  (97) 

/>W 

Hm  =  cpd(Td  -  Tie f)  +  /  (Hv[Tie f]  -  qst[Q,  7’,-el])  d.Q.  (98) 

JO 

The  correlations  were  validated  by  comparing  the  effectiveness  values  obtained  from 
these  correlations  with  effectiveness  values  provided  by  MOSHMX  [43]  and  the  ideal 
dehumidifier  theory  [33].  The  results  showed  that  a  e-NTU  model  incorporating  these 
correlations  can  rapidly  and  accurately  predict  the  dehumidifier  performance.  And  they 
were  suitable  to  be  used  in  long  term  performance  simulation  of  desiccant  cooling  system. 

Also,  by  using  the  same  e-NTU  model,  Klein  et  al.  [1]  analyzed  a  regenerative  enthalpy 
exchanger,  Van  den  Bulck  et  al.  [47]  investigated  a  desiccant  cooling  system  and  a  modified 
method  was  proposed  to  solve  s-NTU  model  for  the  desiccant  which  has  discontinuous 
isotherms  and  isopiestics  such  as  LiCl  [48]. 

4.1.2.  Radial-flow  PGS  model 

The  one-dimensional  desiccant  wheels  focused  on  axial  flow  have  been  extensively 
investigated  in  the  literatures  as  described  in  previous  section.  Elsayed  and  Chamkha  [34] 
studied  the  radial  flow  desiccant  wheel  and  developed  a  model  to  predict  its  performance. 
The  radial  flow  wheel  with  an  inner  radius  rx  and  outer  radius  r2  is  shown  in  Fig.  8.  It  is 


Regeneration 

< 


Dehumidified 

> 


Fig.  8.  Schematic  of  radial-flow  desiccant  wheel  and  the  model  of  control  volume  for  radial-flow  PGS  model  [34]. 
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divided  into  two  sections:  period  1  of  fraction  /?  for  regeneration  air  and  period  2  of 
fraction  1— /?  for  process  air.  The  air  flows  in  the  radial  direction  inwards  or  outwards.  The 
assumptions  for  the  model  were  similar  to  that  of  Jurinak  and  Mitchell  [9]  and  Van  den 
Bulck  et  al.  [46]. 

Conservation  of  moisture  in  the  air  was  expressed  by 


d  AfdYa 
At  dr 


dr  =  KryAv  d At  d r(Ym 


Conservation  of  moisture  in  the  matrix  was  given  by 


2nrdr  0W^ 

mmco—^2 - dt  =  K'vAv  d-4'  Ya  ~  Ym). 

n(j2  —  rf)  dt 

Conservation  of  energy  in  the  air  yielded 
dArdH 


(100) 


m, 


-i 


dr  +  K'Av  d Aj  d rHv(Ym  -  Ya)  =  hAv  d Af  d r(Tm  -  Ta). 


Aj  dr  '  '  y 
Conservation  of  energy  in  the  matrix  yielded 


(101) 


min  o) 


2nrdr  dHm  .  d A{dHg 

dr  =  0. 


7 —  rf)  dt 


dt  +  rh; 


A;  dr 


(102) 


The  boundary  conditions  in  modes  I  and  II  were  applied  in  the  model.  Periodic 
boundary  condition  was  expressed  by  the  following  equations: 


W(r/rucot  =  0)  =  W(r/rucot  =  1),  (103) 

Tm(r/ru  cat  =  0)  =  Tm(r/ruwt  =  1).  (104) 

Based  on  experimental  data,  the  corrected  compound  heat  and  mass  transfer  coefficients  h 
and  K  were  dependent  on  the  superficial  air  velocity  u  within  a  matrix  cross-sectional 
area  [49] 

h'  =  Chua ,  (105) 


K'y  =  Cmuapa,  (106) 

where  C/7,  Cm  and  a  are  constants. 

The  relations  to  determine  the  heat  of  sorption,  equilibrium  humidity  ratio  of  air  in 
contact  with  a  silica  gel  matrix  and  the  enthalpy  of  the  silica  gel  were  the  same  as  Eqs.  (35), 
(96)— (98). 

These  differential  governing  equations  were  reduced  to  dimensionless  NTU  style  and 
discretized  by  finite  difference  technique  on  staggered  grid.  Parametric  studies  were 
performed  using  the  model  to  determine  the  effects  of  flow  configuration  and  operational 
parameters  on  the  performance  of  radial  flow  silica  gel  desiccant  wheel  were  studied.  Also 
axial  flow  desiccant  wheel  was  compared  with  conventional  radial  flow  wheel.  It  was 
founded  that  whose  performance  is  better  depending  on  the  values  of  the  design  and 
operating  parameters.  For  instance,  in  ventilation  and  recirculation  modes,  a  radial  flow 
wheel  with  a  =  2  (Eq.  (105))  has  lower  outlet  humidity  ratio  of  process  air  and  higher 
performance  coefficient  than  axial  flow  wheel.  But  when  a  ^6,  the  results  are  opposite. 

It  should  be  noted  that  composite  transfer  coefficients  h  and  K'y  are  used  to  account  for 
solid  diffusion  resistance  in  PGS  models.  Therefore,  although  second-order  differential 
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terms  are  not  added  in  the  governing  equations  and  the  style  of  the  equation  is  almost  the 
same  as  GSR  model.  The  SSRs  are  indirectly  included.  That  is  why  PGS  model  is  excluded 
from  GSR  model  and  included  in  GSSR  model. 

Since  the  composite  coefficients  are  obtained  from  large  numbers  of  experimental 
results,  it  is  not  convenient  for  studying  a  variety  of  desiccant  materials  and  geometries. 
And  by  now,  only  the  data  for  silica  gel  has  been  obtained.  Hence,  the  reliance  on  such 
excessive  available  data  limits  the  application  of  PGS  model.  Besides,  the  reliability  of  PGS 
model  is  not  good  enough  which  hinders  its  further  development.  Clark  et  al.  [50]  tested  a 
prototype  bed  designed  for  solar  air  conditioning  and  found  the  results  were  in  rather  poor 
agreement  with  predictions  based  on  the  PGS  model,  particularly  after  a  step  change  in 
inlet  air  condition.  They  concluded  that  the  discrepancy  was  due  to  the  shortcomings  of  the 
model,  since  the  SSR  exceeds  the  GSR  under  these  conditions.  The  experimental  result 
reported  by  Clark  et  al.  [50]  was  somewhat  limited  and  imprecise  for  the  bed  with  a  large 
prototype  design.  Hence  Pesaran  [21]  performed  bench  scale  experiments  on  thin  beds  of 
regular  density  silica  gel  with  step  changes  in  inlet  air  humidity.  Their  results  conformed 
that  the  SSR  was  indeed  generally  larger  than  GSR.  Also,  the  results  of  Mei  and  Lavan 
[51]  indicated  that  the  predicted  results  obtained  from  PGS  model  were  acceptable 
compared  with  experiment  results,  but  the  coherence  was  not  very  well.  Furthermore, 
Gheelayagh  and  Gidaspow  [52]  pointed  out  that  the  composite  coefficients  h!  and  K'v  that 
represent  the  SSRs  were  not  constant,  their  value  changed  with  the  time.  So  adopting 
constant  composite  coefficients  cannot  reliably  forecast  the  performance  of  dehumidifier. 

4.2.  Gas  and  solid-side  ( GSS)  model 

Since  PGS  models  adopting  lumped  transfer  coefficients  cannot  provide  very  accurate 
results,  and  there  exists  a  demand  for  suitable  models  containing  SSRs.  With  the  increasing 
knowledge  of  diffusion  mechanism  inside  desiccant,  gas  and  solid-side  (GSS)  models  are 
proposed  in  which  heat  and  mass  transfer  diffusional  terms  in  solid-side  are  directly  added. 
Although  GSS  model  needs  more  time  to  compute  the  additional  second  order  diffusion 
terms  compared  with  PGS  model,  it  is  more  related  to  the  true  physics  of  the  problem. 
Therefore,  GSS  model  is  likely  to  give  better  agreement  with  experiment  results. 

4.2.1.  One- dimensional  model 

Using  the  control  volume  illustrated  in  Fig.  9,  San  and  Hsiau  [35]  have  derived  a  one¬ 
dimensional  model  considering  the  axial  mass  diffusion  and  heat  conduction  in  solid 
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Fig.  9.  The  model  of  control  volume  for  a  one-dimensional  GSS  model  [35]. 
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desiccant,  while  the  storage  terms  of  moisture  and  energy  in  air  dYa/dt  and  dTa/dt  were 
neglected.  The  mass  and  energy  balances  were  expressed  as: 

Conservation  of  moisture  in  air: 


a  Ya  AvV 


m. , 


v 


Ky(Ym  -  Ya). 


a  dz  2  L  y 

Conservation  of  moisture  in  solid: 

a  IT  5  a  2W 

m"4-ff-AvV  2Ds-frT  =  AvVKy{Ya-Y™)- 
where  the  surface  diffusivity  is  evaluated  by 
1.6  x  10"6p, 


Ds  = 


m 
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Conservation  of  energy  in  air: 
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Conservation  of  energy  in  solid: 
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(108) 


(109) 


(110) 


0  T  &  02  rj~i 

mmcpm^-AvV-k^  =  AvVh(Ta  -  Tm)  +  AvVKy(Ya  -  Ym)qst,  (111) 

where  k  is  the  compound  conduction  coefficient  and  calculated  using  the  following 
expression: 


k=f-kd  +  {\  -f)ksp  (112) 

and  subscript  “m”  means  the  solid  part  of  the  wheel  which  contains  desiccant  and 
substrate  materials  and  then  “mm”  stands  for  the  total  mass  of  solid  desiccant  wheel. 

RD  gel  was  selected  as  desiccant  in  the  study.  Due  to  surface  diffusion  is  the  only 
mechanism  need  to  be  considered  for  RD  gel  [39],  only  “Z>s”  which  represents  surface  mass 
diffusion  coefficient  appears  in  the  model.  The  heat  of  sorption  and  equilibrium  isotherm 
relationship  were  determined  by  Eqs.  (35)  and  (96).  The  saturation  pressure  of  vapor  Pvs 
was  evaluated  using  the  following  expressions: 


vs 


218.167 


z  a  +  bz  +  cz 3 

T  1  +  dz 


(113) 


z  =  100.97  —  T;  a  =  3.2437814;  b  =  5.86826  x  10"3, 
c  =  1.1702379  x  10"8;  d  =  2.1878463  x  10"3. 


The  relationship  between  Ym,  Tm  and  <P  was  obtained  from  ASHRAE  data 

6.22  x  10’3<f> 

Y'n  ~  1Q4.21429-[7.5(r,„-273.15)/(3’, ,,-35.85)]  _  Q  ' 


(114) 


A  mixed  numerical  scheme  was  developed  to  solve  the  governing  equations.  The  two 
governing  equations  in  desiccant  side  (Eqs.  (108)  and  (1 1 1))  were  solved  using  full-implicit 
finite  difference  scheme  and  the  other  two  in  gas  side  (Eqs.  (107)  and  (110))  were  solved 
using  backward  finite  difference  scheme.  The  full-implicit  finite-difference  scheme  was 
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verified  to  be  effective  for  solving  the  equation  with  diffusion  terms.  Also  this  mixed 
scheme  was  compared  with  the  explicit  scheme  which  solved  the  four  governing  equations 
totally  with  a  backward  finite  difference  scheme.  It  showed  that  for  the  same  accuracy  the 
mixed  scheme  takes  less  computation  time  than  the  explicit  scheme.  On  the  other  hand,  for 
the  same  time  step  At  the  balance  in  the  mixed  scheme  is  much  better  than  that  in  the 
explicit  scheme. 

Numerical  results  of  this  model  indicated  that  the  axial  heat  and  mass  diffusion  are 
governed  by  two  important  parameters  NTU  and  3  /4 Bi,  and  for  a  case  of  zero  SSR,  the 
calculated  performance  of  the  silica-gel  desiccant  wheel  is  poor.  Furthermore,  it  was  found 
that  only  taking  account  of  surface  diffusion  is  good  enough  for  predicting  RD  silica  gels 
desiccant  wheel,  but  for  ID  silica  gels  the  Knudsen  and  ordinary  diffusion  are  needed  to  be 
considered. 


4.2.2.  Two-dimensional  model 

Charoensupaya  and  Worek  [10]  established  a  two-dimensional  GSS  model  in  which  the 
heat  conduction  and  mass  diffusion  including  gas  phase  diffusion  and  surface  diffusion  in 
the  radial  direction  within  solid  desiccant  were  considered.  The  control  volume  is 
illustrated  in  Fig.  10. 

Conservation  of  moisture  in  the  process  air  was  expressed  as 


1  QYg  dYg\ 
u  dt  dz  ) 


=  2  Ky(Yd 


(115) 


Conservation  of  moisture  in  the  desiccant  felt  was  given  by 


dYd  ,,  ,  dW  n  d2Yd  n  d2W  n 

^ tPad  rs.  f  0  &t)Pd  ~EZ  ^ GP ad  Pd  ^  o  0, 


dt 


dt 


dr2 


dr2 


(116) 


where  DG  and  Ds  represent  the  effective  gas  phase  diffusivity  and  surface  diffusivity, 
respectively. 

Conservation  of  energy  within  the  desiccant  felt  gave 


SHad 

Z'Pad 


+  (1  —  £t)PdCpd 


87V 

dt 


=  (1  -  £t)pdqs, 


8  W 

OF’ 


(117) 


where  the  subscript  “ad”  means  the  air  in  desiccant  pore. 


nonhyroscopic  wall 


L _ I  - ►  Z 


Fig.  10.  The  model  of  control  volume  for  a  two-dimensional  GSS  model  [10]. 
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The  rate  of  energy  transfer  between  process  air  and  desiccant  felt  yielded 


dH 


9 


dt 


+  u 


d  Hn  2  K 


9 


y 


dz 


PgU 


(Yd(t,X,  Ctd)  -  Y a) 


dH 


9 


a 


a  y 


+ 


2  h 


a 


PgU 


(Td(t,x,  otd)  -  Ta). 


(118) 


a 


The  boundary  conditions  were  in  the  form  of  mode  III.  Also,  the  temperature  of  the  wall 
was  assumed  to  be  the  same  as  the  desiccant  temperature  at  r  =  0 


Tw(t,z)  =  Td(t,z,  0).  (119) 

The  model  was  verified  using  the  following  equilibrium  relation  and  heat  of  sorption: 

W  =  0.50,  (120) 


qst  =  hfg  (121) 

The  heat  and  mass  transfer  coefficients  were  approximately  by  constant  value  that  can  be 
calculated  from  the  Nusselt  and  Sherwood  number  correlations. 

An  explicit  finite-difference  scheme  was  adopted  as  the  solution  method  to  calculate  the 
model.  The  model  was  used  to  analyze  the  effect  of  some  parameters,  including  heat  and 
mass  transfer  Biot  number,  on  the  performance  of  desiccant  cooling  system.  It  was 
revealed  that  there  exist  optimum  values  of  COP  and  cooling  capacity  for  given  operating 
conditions.  The  optimum  values  can  be  reached  by  controlling  the  system  cycle  time.  Also 
the  accuracy  of  this  model  has  been  verified  by  experimental  data  in  Ref.  [53]. 

Yu  et  al.  [36,54]  also  constructed  a  two-dimensional  mathematical  model  for  desiccant 
wheel  based  on  the  control  volume  shown  in  Fig.  2.  Both  the  mass  diffusion  and  heat 
conduction  inside  the  solid  desiccant  wheel  in  axial  as  well  as  circumferential  directions 
were  considered  in  the  model.  The  governing  equations  are  listed  as  follows: 

Conservation  of  moisture  in  the  air : 


dYg  dYg  dYg 

dt  dcp  pjFs  dz 


KyAv 


(Yd 


Conservation  of  moisture  in  solid  adsorbent : 


(122) 


8W 

dt 


“I-  co 


dW 

dcp 


2  ln(r2/fi)02  W  3  2W' 

r2  ~  r\  9z2 


KyAv 

Pd 


( Ya 


(123) 


Conservation  of  energy  in  the  air : 


dTg  dTg  m^_dTg 

dt  dcp  PjFs  dz 


PiFsippg  T 


YnC 


a^pv 


) 


Conservation  of  energy  in  solid  desiccant : 


(124) 


dTd  dTd  kd{\  —  Fs)  \2\n(r2/rx)d2Td  d2Td 

dt  dcp  [pd(cpd  +  Wcpi)  +  [  r22-r\  dcp2  dz 2 

_  [hAATg  -  Tf>  +  KyAy(Y a  -  YpqJ 

\-Pddpd  4“  Wc,,i )  +  PspCpsp\ 


(125) 


The  boundary  and  initial  conditions  were  similar  to  those  of  Ref.  [14].  Eq.  (26)  was 
adopted  to  calculate  Yd.  Equilibrium  isotherms  0,  obtained  by  fitting  fourth-degree 
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polynomials  to  the  manufacture  data  [21],  and  heat  of  sorption  qst  for  RD  and  ID  were: 
For  RD  gel : 

<P  =  0.0078  -  0.05759  fF  +  24.16554  W2  -  124.478  W3  +  204.226  W4,  (36) 


—  12400 IF  +  3500  1FsS0.05  (kJ/kg  water), 

—  1400  IF  +  2950  IF  >0.05  (kJ/kg  water), 


For  ID  gel : 


1.235 IF  +  267.99 IF2  -  3170.7JF3  +  10087.161F4  Ws;0.07, 
0.3316  +  3.181F  W>0.07, 


(126) 


(127) 


-300W  +  2095  Wsc0.15  (kJ/kg water), 
2050  W>0.15  (kJ/kg water). 


(128) 


Central  difference  scheme  was  applied  to  discretize  the  diffusion  items  and  the  same 
schemes  as  in  Ref.  [8]  were  adopted  to  discretize  the  other  terms.  It  was  demonstrated  the 
numerical  results  of  this  model  are  very  close  to  the  steady  and  transient  experimental  data 
of  Refs.  [55,56]. 

Dai  et  al.  [57]  pointed  out  that  the  model  proposed  by  Yu  et  al.  [36]  was  helpful  for 
predicting  the  performance  of  desiccant  wheel  but  was  not  brief  enough.  Hence  they 
adopted  psychrometric  wave  chart  analysis  associated  with  the  finite  difference  model 
trying  to  obtain  the  numerical  results  efficiently.  In  their  study,  RD  silica  gel  was  used  as 
adsorbent.  According  to  the  conclusion  of  Pesaran  [39],  the  surface  diffusion  is  dominant 
for  RD  silica  gel,  and  then  the  effective  diffusivity  was  given  as 

-0.974  x  10"6  x  E- ) ,  D0  =  0.8  x  10“6.  (129) 

Using  the  wave  analysis  in  psychrometric  chart  incorporating  with  numerical 
calculation,  the  influences  of  some  important  parameters,  such  as  heat  capacity, 
adsorption  heat,  rotation  speed,  thickness  of  the  desiccant  matrix,  etc.,  on  the  performance 
of  RD  silica  gel  desiccant  wheel  system  were  discussed.  It  was  found  that  the  chart  method 
is  rapid  in  evaluating  the  system  performance  and  is  useful  in  system  design. 

Gao  et  al.  [58]  have  described  a  mathematical  model  that  can  be  seen  as  the 
simplification  form  proposed  by  Yu  [36],  where  diffusion  terms  in  circumferential  direction 
were  ignored  and  the  sensible  heat  transfer  was  included  in  the  governing  equations.  They 
conducted  a  parameter  analysis  employing  the  model  to  investigate  the  effects  of  passage 
shape  and  geometrical  size  of  the  matrix  on  the  performance  of  desiccant  wheel. 

Zhang  and  Niu  [37]  have  presented  a  two-dimensional  transient  model  in  which  heat 
conduction,  surface  and  gaseous  diffusion  in  axial  as  well  as  radial  directions  of  the 
desiccant  wheel  are  included  simultaneously.  Due  to  symmetry,  the  mid-plane  of  a  channel 
was  treated  as  adiabatic  wall,  and  a  half-size  channel  was  chosen  as  the  control  volume  as 
shown  in  Fig.  11.  The  derived  governing  equations  were: 

Conservation  of  moisture  in  the  air  stream : 


De  =  Ds  =  D0  exp 


ldYa  QYa_  4 Ky 

+  T —  —  7/ - ( 

u  ot  oz  aeupg 


(130) 
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Fig.  11.  The  model  of  control  volume  for  a  two-dimensional  GSS  model  [37]. 


Conservation  of  moisture  in  the  desiccant : 


^tPad 


QYd 
dt 


+  Pd 


a  w 
a  t 


Pad 


dz 


+  Pd 


a  /  a  ya  a  /  a  ya 
G~di~ )  +a^ (  ). 

a  /  a w\  a  /  sir 
az  V  s  az )  ar  V  s  a r 


(131) 


where  the  effective  surface  diffusivity  Ds  can  be  calculated  by  Eq.  (129)  and  DG  which 
stands  for  the  composite  effective  diffusivity  of  Kundsen  and  ordinary  diffusion  is  given  by 


Dg  =  ~ 


1  1 
+ 


D0  =  1.758  x  10 


7-1.685 
-4  1  d 

B 


Z>k  =  97  a 


(-) 

\yyi\J 


0.5 


(132) 


T  \Dq  Z>k 

where  a  (m)  is  the  pore  radius  of  the  adsorbent  and  m\  (kg/mol)  is  the  molecule  weight  of 
water. 

Conservation  of  energy  in  the  air  stream : 


\dTa  +  dTa 


4  h 


(Ts  -  Ta). 


u  dt  a z  deupgcpg 
Conservation  of  energy  in  the  desiccant : 

dTd  .  /0 2Td  d2TA  dW 


(133) 


Pdcpt 


a  t 


,  (CTd  CTd\ 

~kd\5E  +  EE)  +qs,Pd 


dt 


(134) 


where  cpt  is  the  total  heat  capacity  of  wet  desiccant,  which  included  two  parts:  dry 
desiccant  and  adsorbed  water  and  is  calculated  by 

Cpt  =  Cpd  +  Wcpj.  (135) 

By  using  heat  mass  transfer  analogy,  the  relations  between  h  and  Ky  can  be  expressed  as 

Ky  =  -f-  ■  (136) 

CpgEe 

Water  content  in  the  desiccant  is  governed  by  a  general  sorption  isotherm  as 

fw . 


w  = 


max 


1  -R  +  R/& 


(137) 


The  boundary  conditions  were  similar  to  Eqs.  (5)— (8),  (12),  (14).  This  model  was 
normalized  and  numerically  solved  by  means  of  alternating  direction  implicit  (ADI) 
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method  [59].  Using  the  model,  effects  of  rotary  speed,  NTU  and  specific  area  on  the 
performance  of  silica  gel  wheel  were  studied  and  compared  in  desiccant  wheels  for 
dehumidification  and  enthalpy  recovery. 

Niu  and  Zhang  [60]  performed  an  analysis  study  adopting  the  above  model  to  study  the 
influence  of  wall  thickness  on  the  optimum  rotary  speeds  for  dehumidification  desiccant 
wheel  and  enthalpy  recovery  wheel. 

Sphaier  and  Worek  [38]  derived  a  desiccant  wheel  mathematical  model  which 
accounts  for  heat  and  mass  transfer  processes  inside  solid  desiccant  too.  The  model  is 
different  from  others  in  two  aspects:  (1)  the  diffusion  terms  in  the  governing  equations 
are  given  in  the  form  of  vector,  hence  governing  equations  can  be  easily  converted 
into  one  or  two,  even  three-dimensional  ones;  (2)  parameters  describing  detailed 
constitution  of  felt  are  included  in  the  model.  The  model  is  thought  to  be  useful  for 
analyzing  the  transport  phenomena  within  the  porous  adsorbent  felt  and  the  influence  of 
supporting  structure.  It  is  assumed  that  the  entire  felt  (subscript  “/’)  is  composed  of  rigid 
matrix  (subscript  “fm”)  and  pores  (subscript  “fp”)  in  which  both  the  gas  and  adsorbed 
liquid  phase  coexist.  The  matrix  includes  solid  adsorbent  desiccant  (subscript  “/if’)  and 
inert  materials  which  has  no  effects  on  the  sorption  characteristics  of  adsorbent  and 
surface  diffusion. 

Based  on  the  control  volume  shown  in  Fig.  12,  the  derived  conservation  equations  of 
moisture  and  energy  in  the  air  stream  were  similar  to  Eqs.  (130),  (133),  and  conservation 
equations  of  moisture  and  energy  in  the  desiccant  were: 

dYf  0  Wf 

£,pa  “0T  +  (1  “  e,)V-0T  =  PaV(DG,eKSYfp)  +  Pyi/"V(Z)Sjeff V  Wfm),  (138) 


dT 


f 

PfCf  ~kt  =  V(V,effV7»  +  pfdf 


dt 


(1  ~St)~dfL~  v(DSfiaswfm) 


tfsf 


(139) 


where  the  effective  heat  conductivity  kf  Q ff  gas  phase  and  surface  diffusivity  DG  q ff,  DS  e ff  are 
defined  as 


^7,e ff  —  Zfkfp  T  (1  £ t)kfm ? 


(140) 


Dg,q  ff  =  (141) 

A>,eff  =  (1  -  St)DS/Ts.  (142) 

Eqs.  (7)  and  (8)  were  utilized  to  depict  the  insulated  and  impermeable  boundaries. 
Boundary  conditions  in  mode  III  were  obtained  through  interface  balances:  at  r  =  Rg 
(Fig.  12),  the  boundary  condition  for  mass  transport  was  the  same  as  Eq.  (12)  and  for  heat 


Impermeable  wall,  w 

-Z -  r=R 


Sorbent  felt,  f 

- r 

i 


r=R 


i 


r=R 
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Process  stream,  g 


I  Az 

,< — > 


Fig.  12.  The  model  of  control  volume  for  a  GSS  model  with  vector  form  [38]. 


1520 


T.S.  Ge  et  at  /  Renewable  and  Sustainable  Energy  Reviews  12  (2008)  1485-1528 


transport  was  given  as: 

—  ^/,eff  =  h(T a  —  Tf)  +  Ky(Y  a  —  Yfp)q^T  —  Py#7)s,eff  0^  list 

and  at  the  wall  the  boundary  conditions  were 


(143) 


a  y 


~PaDG,e ff  ~  PfdD^fi  — ^ 


3  W 


fin 


=  0. 


(144) 


•  _  /  n  ®  Yfp  dWfin 

qs  ~  -*U ff_07  “  PaDG,eff~^rqAT  ~  P/^S,eff  H  / , 

where  qAT  is  termed  the  heat  of  sorbate  transfer  at  the  interface 


(145) 


qAT  =  -  ^).  (146) 

Also  the  selecting  equilibrium  isotherm  was  in  the  form  of  Eq.  (17). 

These  governing  equations  were  fully  normalized  using  classical  dimensionless  groups 
for  heat  and  mass  transfer,  and  a  numerical  algorithm  based  upon  the  finite  volumes 
method  was  employed  to  solve  it.  It  showed  that  simulation  results  present  a  very 
reasonable  agreement  with  the  experimental  data.  Moreover  the  simulation  results  of  a 
test-case  were  presented,  suggesting  a  possible  optimization  to  wheel  initial-cost  and 
compactness  by  reducing  felt  thickness.  This  model  can  be  applied  for  the  simulation  of 
several  configurations  such  as  desiccant  wheels  or  enthalpy  exchangers. 


4.3.  Parabolic  concentration  profile  (PCP)  model 

Pseudo-gas-side  (PGS)  model  and  gas  and  solid-side  (GSS)  model  have  their  own 
shortcomings.  The  GSS  model  requires  greater  computational  effort  than  PGS  model 
because  of  the  additional  second-order  heat  and  mass  transfer  diffusional  items.  However, 
the  PGS  model  requires  extensive  experimental  data  for  determination  of  the  lumped 
transfer  coefficients  with  different  desiccant  materials  and  its  reliability  is  not  good 
enough.  Hence,  a  computationally  efficient  model  that  effectively  considers  both  the  SSRs 
and  GSRs  in  desiccant  wheel  is  needed.  Parabolic  concentration  profile  (PCP)  model 
assuming  an  analytically  convection  concentration  profile  within  the  desiccant  particle  to 
account  for  solid-side  diffusion  is  developed  to  achieve  this  objective.  Usually,  based  on 
the  assumed  concentration  profile,  a  differential  governing  equation  describing  mass 


Fig.  13.  Schematic  figure  of:  (a)  Diffusion  through  a  spherical  particle,  (b)  A  side  view  of  one  of  ducts  for  a  PCP 
model  [39]. 
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diffusion  within  a  desiccant  particle  is  developed  and  included  in  PCP  model  to  represent 
the  moisture  conservation  in  solid  desiccant. 

Pesaran  and  Mills  [39,61]  proposed  a  PCP  model  for  simultaneous  heat  and  mass 
transfer  in  a  thin  packed  bed  of  silica  gel  particles.  It  believed  that  diffusion  of  moisture 
through  the  porous  particle  consists  of  Knudsen  diffusion  and  surface  diffusion.  The 
spherical  particle  model  and  adsorption  bed  model  are  illustrated  in  Fig.  13. 

The  general  form  of  moisture  conservation  equation  for  the  spherical  particle  with  radial 
symmetry  was  derived  and  given  by: 


dW  _  1  6 
dt  r2  dr 


(147) 


where  De  is  the  total  effective  diffusivity  which  includes  Knudsen  diffusion  and  surface 
diffusion: 


a'(W )  fdY 

De  =  £*S,eff  +  Ac.eff  - - ,  d'i  W)  =  Pg  ( 


(148) 


For  the  thin  beds,  storage  term  dYg/dt  was  negligible  and  moisture  conservation  equation 
in  the  air  was 


0T 


m, 


g 


=  Ky(Ys  -  Yq\  1  -  Yq)P . 


vg  0Z  —  s  1  gj\ A  '  g , 

Conservation  of  energy  in  gas-side  was  described  by 

0  T, 


(149) 


cpgmg 


g 


dz 


=  -P[h  +  .  Ky(YS  -  Yg)]  .  (Tg  ~  T S) 


(150) 


where  axial  and  radial  conduction  and  the  storage  term  dTg/dz  have  been  neglected  and 
the  bed  is  assumed  to  be  adiabatic. 

Energy  conservation  in  the  solid  desiccant  yielded: 


APjwC 


rn^pm 


dl\ 

dt 


=  P[h(Tg  -  Ts)  -  qstKy(Ys  -  Yg)]9 


(151) 


where  subscript  “m”  represents  the  adsorption  bed. 

Eqs.  (147)  and  (149)  were  linked  by  the  equilibrium  relation  at  the  particle  surface 


Ys(z,  t)  =  [W{r  =  R,z,  0,  Ts(z,  t),Ps]. 

The  boundary  conditions  and  initial  condition  for  Eq.  (147)  were 


(152) 


dW 


dr 


=  0. 


(153) 


r= 0 


dW 


Pparticle-^  g^,  r=R  —  V  [  Ys(z,  t)  Y g{z,  t)\. 


(154) 


W(z,r,t  =  0)=  W0.  (155) 

The  boundary  condition  for  Eq.  (150)  was 
Tg(z  =  0  ,t)=  Tm 


(156) 
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While  the  initial  condition  for  Eq.  (151)  was 

Ts(z,t  =  0)=T0.  (157) 

Heat  and  mass  transfer  coefficients  were  obtained  by  experimental  correlations  [21]: 

Ky  =  1.70m^Re“°-42,  (158) 

h=  l.60maRe-°42cpg.  (159) 

The  properties  for  RD  gel  and  ID  gel  were  the  same  as  Eqs.  (36),  (1 26)— (128),  and  specific 
heats  of  the  gas  and  felt  were 

cpg  =  1884  Yg  +  1004(1  -  Yg\  (160) 


r  R 

cpm  =  4186PFavg  +  921  Wan=  4nr2WpP  dr/(4/3)nR3 pP.  (161) 

Jo 

The  governing  equations  were  put  into  dimensionless  form  and  solved  by  finite 
difference  method.  Crank-Nicholson  method  was  adopted  for  Eq.  (147).  Implicit  Euler 
method  was  applied  for  Eq.  (151).  A  fourth-order  Runge-Kutta  technique  was  used  for 
Eqs.  (149),  (150).  The  transient  results  of  this  model  were  compared  with  the  predictions  of 
PGS  model.  It  was  proved  that  the  simulated  results  getting  from  the  new  model  agree 
better  with  the  experimental  data  than  those  from  PGS  model  for  both  ID  gel  and  RD  gel 
desiccant  material. 

Chant  and  Jeter  [40]  performed  a  steady  state  and  transient  analysis  of  rotary  desiccant 
dehumidifier  using  a  PCP  model.  The  heat  transfer  equation  and  conservation  equation  in 
the  desiccant  side  were  similar  to  those  of  PGS  model  (Eqs.  (77),  (78),  (81)),  except  for  gas- 
side  transfer  coefficients  h  and  Kv  not  compound  transfer  coefficients  like  h  and  K'v  were 
adopted  in  PCP  model.  For  mass  transfer  analysis,  a  parabolic  concentration  profile  was 
assumed  to  exist  in  desiccant  particle  at  all  times,  and  for  a  spherical  particle  the  profile 
was 


W  =  ci\  +  ci2(r/R)2. 

The  moisture  conservation  equation  for  a  spherical  particle  was 


dW 

dt 


3  Ky 


P  particle^ 


[Yd(Ws,  Td) 


(162) 

(163) 


where  Ws  is  desiccant  surface  concentration  which  is  derived  using  Eq.  (162),  Yd  stands  for 
moist  air  equilibrium  humidity  ratio,  they  could  be  calculated  by  the  following  equation: 


Ky  ^particle 
^Pparticle^e 


[Yd(Ws,  Td) 


Ya]  +  W, 


(164) 


Yd  —  Co  +  c\  W  +  C2Hd  +  C3  Y a.  (165) 

[62]  The  enthalpy  and  equilibrium  adsorption  isotherm  expressions  were  similar  to 
Eqs.  (97)-(98)  and  Eq.  (17). 

First  the  governing  equations  (Eqs.  (77)-(78),  (81),  (163))  were  associated  with  the  inlet 
states  of  moist  air  for  processing  and  regeneration  as  well  as  the  desiccant  bed  initial  state 
to  construct  transient  model,  which  was  used  to  study  startup  effects.  Then  periodic  steady 
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state  model  was  distinguished  by  applying  the  following  periodic  boundary  conditions: 


Hd(z,  t  =  0)  =  Hd(z,  t  =  tr). 

(166) 

5 

FT 

II 

o 

II 

FT 

II 

(167) 

Transient  model  and  periodic  steady  state  model  constituted  the  whole  PCP  model.  The 
transient  model  was  solved  using  Bulirsch-Stoer  integrator  method  [63].  An  implicit 
periodic  steady-state  solution  with  a  linear  adsorption  isotherm  (Eq.  (165))  which  utilizes  a 
sparse  matrix  solving  routine  was  developed  to  solve  periodic  steady-state  model. 

Compared  with  the  transient  experimental  results  of  a  fixed  bed  dehumidifier.  It  was 
found  that  the  PCP  model  does  not  predict  well  in  the  initial  adsorption  period,  which  is 
due  to  the  reason  that  the  moisture  concentration  profile  is  not  parabolic  until  later  in  the 
sorption  period.  But  the  simulation  results  become  similar  to  that  of  the  PGS  model  as 
time  moves  on.  Also  the  periodic  steady  state  of  a  desiccant  wheel  system  predicted  by  this 
numerical  model  has  shown  good  agreement  with  the  testing  data  of  practical  desiccant 
wheel  at  low  and  medium  rotation  rates  condition.  The  authors  discussed  the  efficiency 
and  flexibility  of  PCP  model  too.  Their  results  indicated  that  the  model  only  requires 
1-3  min  on  a  25  MHz  80486DX  PC  for  the  cased  considered.  And  since  the  widely  used 
adsorbent  property  effective  diffusivity  De  is  used  in  mass  transfer  equation,  it  has  the 
capability  to  simulate  a  variety  of  desiccant  materials.  Hence,  PCP  model  is  extremely 
efficient  and  an  excellent  tool  for  investigating  alternative  desiccant  materials. 

5.  Empirical  model 

Beccali  et  al.  [64]  believed  the  models  based  on  the  detailed  heat  and  mass  transfer 
phenomena  inside  desiccant  wheel  are  too  complex.  The  solving  process  is  complicated  and 
time-consuming.  Hence,  in  order  to  evaluate  the  performance  of  various  kinds  of  rotary 
desiccant  wheels  available  in  the  market,  they  analyzed  the  available  measurements  data 
from  the  manufacturer  thoroughly  and  developed  empirical  model  for  evaluating  the 
performance  of  silica  gel  desiccant  wheel.  It  should  be  emphasized  that  the  units  of 
enthalpy,  temperature  and  humidity  ratio  are  KJ/kg,  °C  and  g/kg,  respectively,  in  the 
study.  The  derived  correlations  were  expressed  as 

X  out  /  T  out  —  (Fx  Yl  +  F2  Ym  +  F3)Tm  +  (F4  Yl  +  F5  Yin  +  F6),  (168) 

where  Ft  is  a  second-order  polynomial  of  Freg  and  can  be  expressed  as  follows: 

i  =  {AjT2YQg  +  BjTveg  +  C/),  i  =  1, . . . ,  6,  (169) 

where  Ah  Bh  Q  are  second-order  polynomials  of  Xreg, 

Aj,  Bj,  Cj  =  at  F^eg  +  bj  Y  reg  +  Cj,  i=l,.  ..,6.  (170) 

The  equations  were  used  to  forecast  process  air  outlet  temperature  rout  and  absolute 
humidity  Fout  as  a  function  of  the  temperature  and  humidity  ratio  of  the  regeneration  air 
(Freg  Freg)  and  of  the  inlet  state  of  process  air  (T-m  Fin).  There  are  54  parameters  (a,  bf  ct) 
that  related  to  outlet  condition  Fout  and  Fout  in  the  model,  so  the  model  is  named  as 
“Model  54”.  Model  54  was  carried  out  to  predict  the  performance  of  silica  gel  desiccant 
wheel  (type  -I).  It  was  observed  that  the  calculated  values  obtained  using  this  model  are 
very  close  to  the  measured  values. 
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However,  for  a  different  type  of  desiccant  material,  a  new  group  of  54  coefficients  are 
needed,  which  is  a  tedious  job  for  users.  Then,  further  simplification  were  attempted  to 
overcome  the  complexities  of  “Model  54”.  The  new  model  was  expressed  by  the 
parameters  0  (relative  humidity)  and  H  (enthalpy) 


A0  —  (^in  0out)  —  Cl(0 in  ^reg)  T  b , 

(171) 

AH  =  (Hm  -  Houl)  =  a'(Hreg  -  Hm)  +  b', 

(172) 

where  the  enthalpy  is  calculated  by 

(2501  +  1.85  T)Y 

H=  .  7  +  1.0067 

1000 

(173) 

the  calculation  of  relative  humidity  ratio  is  based 
Beccali  [65] 

on  empirical  relation  developed  by 

<P  =  (18.67157+  1.7976)e-°-0537’. 

(174) 

The  new  relational  expressions  only  have  four  variables  a ,  b ,  a ',  b\  much  less  than  the 
variable  numbers  of  “Model  54”.  Three  different  kinds  of  desiccant  wheels  were  used  to 
validate  this  model,  two  are  based  on  silica  gel  and  another  one  is  on  lithium  chloride.  The 
results  showed  that  the  simplified  model  reasonably  predict  the  performance  of  all  wheels. 

The  above  two  empirical  models  are  valid  only  for  the  wheel  running  with  identical  air 
flux  in  process  and  regeneration  side.  By  adding  correction  factors,  Beccali  [66]  further 
modified  the  model  which  can  be  used  to  forecast  the  system  with  unequal  process  and 
regeneration  air-flow  volume.  The  modified  model  was  verified  for  the  three  different  kinds 
of  desiccant  wheels. 

6.  Conclusion 

Various  mathematical  models  have  been  proposed  to  assess  the  performance  of 
desiccant  wheel  for  different  operating  design  conditions.  Generally,  the  mathematical 
models  are  based  on  heat  and  mass  transfer  phenomena  inside  the  system.  They  can  be 
divided  into  two  main  categories  in  terms  of  gas-side  resistance  (GSR)  model  and  gas  and 
solid-side  resistance  (GSSR)  model.  Within  the  two  categories,  the  models  can  be 
subdivided  based  on  the  dimensions  and  the  methods  to  represent  solid-side  resistance 
(SSR).  GSR  model  has  two  dimension  types:  one-dimensional  and  two-dimensional. 
GSSR  model  can  be  classified  into  three  sub-categories:  pseudo  gas-side  (PGS)  model,  gas 
and  solid-side  (GSS)  model  and  parabolic  concentration  profile  (PCP)  model.  Also,  using 
the  measurements  data,  several  empirical  models  to  evaluate  the  performance  of  desiccant 
wheel  are  reviewed  in  this  paper. 

In  the  GSR  model,  heat  and  mass  transfer  within  solid  desiccant  is  not  taken  into 
account.  The  governing  equations  have  relatively  simple  form  with  lower  accuracy.  PGS 
model  adopts  the  lumped  heat  and  mass  transfer  coefficients  to  represent  the  overall  heat 
and  mass  transfer  from  the  air  stream  to  desiccant  material.  In  theses  models,  convective 
heat  and  mass  transfer  coefficients  are  empirically  degraded  to  account  for  both  gas-side 
and  solid-side  resistances.  Therefore,  though  its  form  is  basically  the  same  as  GSR  model, 
the  accuracy  is  improved.  The  more  practical  model  known  as  GSS  model  adds  second- 
order  heat  and  mass  transfer  diffusional  terms  directly  into  the  governing  equations.  It  is 
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more  related  to  the  actual  process  occurred  in  the  desiccant  wheel  and  the  effect  of  SSRs 
can  be  well  explained.  The  precision  of  this  model  is  thus  greatly  improved.  PGS  model 
consumes  less  computational  time  but  has  lower  flexibility  and  accuracy  compared  to  GSS 
model.  PCP  model  which  assumes  an  analytically  convection  concentration  profile  within 
the  desiccant  particle  appears  to  be  a  good  compromise  between  fundamental  principles, 
accuracy  and  computational  speed.  It  is  an  excellent  tool  for  performing  seasonal 
simulations  or  repetitive  design  simulation  particularly  when  alternative  materials  are 
investigated. 

It  can  be  seen  that  boundary  conditions  in  Mode  I  are  widely  used  in  all  the  models 
whenever  possible.  Mode  II  and  III  boundary  conditions  are  usually  utilized  in  GSSR 
model.  Other  auxiliary  relations  such  as  heat  of  sorption  and  equilibrium  isothermal  have 
to  be  obtained  to  complete  the  mathematical  model. 

Due  to  the  complex  nature  of  the  partial  differential  equations,  it  is  necessary  to  solve 
the  mathematical  model  using  numerical  methods.  Finite  difference  method  is  recognized 
as  the  most  accurate  and  most  universal  solution  technique.  Though  other  methods  such  as 
finite  volume  method  [7,38],  analogy  method  [31]  and  wave  analysis  method  [33,57]  are 
also  used. 

Beside,  model  validation  is  an  important  step  in  model  development  since  it  offers  the 
possibility  of  comparing  simulation  results  with  actual  system  behavior.  Experiments  are 
mostly  used  to  validate  the  mathematical  model.  Beside,  compared  with  the  previous 
numerical  results  is  also  a  good  validation  method. 

Though  a  large  amount  of  works  have  been  conducted  on  modeling  and  analyzing 
desiccant  wheel,  further  efforts  are  still  needed: 

(i)  To  take  account  of  the  pressure  and  velocity  loss  through  the  channel. 

(ii)  To  study  the  influence  of  variable  thermophysical  physical  properties  of  moist  air  and 
desiccant  materials  which  are  taken  as  constant  in  almost  all  existing  models. 

(iii)  To  improve  the  accuracy  of  the  model  based  on  diffusion  mechanism,  since  by  now, 
there  still  no  uniform  theory  exists. 

(iv)  To  construct  a  simulation  package  of  the  whole  desiccant  air-conditioning  system  by 
combining  the  model  of  the  wheel  and  other  components  in  the  system.  Using  object- 
oriented  languages,  Chau  and  Worek  [67]  simulated  an  open-cycle  desiccant  cooling 
system  operating  in  the  recirculation  mode.  As  presented  in  the  paper,  the  simulation 
package  is  ongoing  and  in  the  near  future  it  is  hoped  that  a  graphical  user  interface  can 
be  applied  in  it  (Tables  2  and  3). 
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